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Introductory Remarks. 


The former Essay* contained a general method for reducing all the most important pro- 
blems of dynamics to the study of one characteristic function, one central or radical relation. 
It was remarked at the close of that Essay, that many eliminations required by this method in 
its first conception might be avoided by a general transformation, introducing the time ex- 
plicitly into a part S of the whole characteristic function V; and it is now proposed to fix the 
attention chiefly on this part S, and to call it the Principal Function. The properties of this part 
or function S, which were noticed briefly in the former Essay,} are now more fully set forth; 
and especially its uses in questions of perturbation, in which it dispenses with many laborious 
and circuitous processes, and enables us to express accurately the disturbed configuration of 
a system by the rules of undisturbed motion, if only the initial components of velocities be 
changed in a suitable manner. Another manner of extending rigorously to disturbed motion the 
rules of undisturbed, by the gradual variation of elements, in number double the number of the 
coordinates or other marks of position of the system, which was firstinvented by Lagrange, t and 
was afterwards improved by Poisson, is considered in this Second Essay under a form perhaps 
a little more general; and the general method of calculation which has already been applied to 
other analogous questions in optics and in dynamics by the author of the present Essay, is now 
applied to the integration of the equations which determine these elements. This general method 
is founded chiefly on a combination of the principles of variations with those of partial differ- 
entials, and may furnish, when it shall be matured by the labours of other analysts, a separate 
branch of algebra, which may be called perhaps the Calculus of Principal Functions;§ because, 
in all the chief applications of algebra to physics, and in a very extensive class of purely mathe- 
matical questions, it reduces the determination of many mutually connected functions to the 
search and study of one principal or central relation. When applied to the integration of the 
equations of varying elements, it suggests, as is now shown, the consideration of a certain 
Function of Elements, which may be variously chosen, and may either be rigorously determined, 
or at least approached to, with an indefinite accuracy, by a corollary of the general method. 
And to illustrate all these new general processes, but especially those which are connected with 
problems of perturbation, they are applied in this Essay to a very simple example, suggested 
by the motions of projectiles, the parabolic path being treated as the undisturbed. As a more 
important example, the problem of determining the motions of a ternary or multiple system, 
with any laws of attraction or repulsion, and with one predominant mass, which was touched 
upon in the former Essay, is here resumed in a new way, by forming and integrating the 
differential equations of a new set of varying elements, entirely distinct in theory (though little 
differing in practice) from the elements conceived by Lagrange,||and having this advantage, that 
the differentials of all the new elements for both the disturbed and disturbing masses may be 
expressed by the coefficients of one disturbing function. 


* [P. 103.] + [Pp. 160, 161.] 

į (Lagrange, “Sur la théorie générale de la variation des constantes arbitraires,” Mémoires de l’ Institut (1808), 
p. 257; Méc. Anal. 3rd ed. Tome 1, pp. 299-320; Poisson, “Sur la variation des constantes arbitraires dans les 
questions de mécanique,” Journal de l’ Ecole Polyt. Tome vm (1809), p. 266.] 

§ [For the development of the Calculus of Principal Functions see pp. 297-410 of this volume. A brief account 
of the theory is given on pp. 408-410.] 

|| [For Lagrange’s elements see Méc. Anal. 3rd ed. Tome u, Chap. 1, § 2.) 
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Transformations of the Differential Equations of Motion of an Attracting or 
Repelling System. 


1. It is well known to mathematicians, that the differential equations of motion of any 
system of free points, attracting or repelling one another according to any functions of their 
distances, and not disturbed by any foreign force, may be comprised in the following formula: 

X.m(x"da+y" dy+2"5z)=8U: (1.) 
the sign of summation È extending to all the points of the system; m being, for any one such 
point, the constant called its mass, and x, y, z being its rectangular coordinates; while x”, y”, z” 
are the accelerations, or second differential coefficients taken with respect to the time, and 
ôx, dy, dz are any arbitrary infinitesimal variations of those coordinates, and U is a certain 
force-function, introduced into dynamics by Lagrange, and involving the masses and mutual 
distances of the several points of the system. If the number of those points be n, the formula (1.) 
may be decomposed into 3n ordinary differential equations of the second order, between the 


coordinates and the time, 


Loo BOs òU òU 
Mets = i MY = sy mM; eek : (2.) 


and to integrate these differential equations of motion of an attracting or repelling system, or 
some transformations of these, is the chief and perhaps ultimately the only problem of mathe- 
matical dynamics. 


2. To facilitate and generalize the solution of this problem, it is useful to express previously 
the 3n rectangular coordinates x, y, z as functions of 3n other and more general marks of position 
N1» N2» +++ Nan; and then the differential eq uations of motion take this more general form, discovered 
by Lagrange,* 


da? 37 _8U ” 
dt dn; Sn; 8,’ i 

in which 
T =42.m (x? +y’? -+z?). (4.) 


For, from the equations (2.) or (1.), 
òU ( DOE MOYI = 
ce Ss A E cre he lc ae 
òn; d òn; òn; 
(5.) 


> 


5 m (e pat Feet ge) 
rok ae Oe oat S 
in which 


aa 


Sy al ST (6.) 
7+2 CEON hoe 


an’ 
=3.m(2'554 $ ; 
Y ni Sn; 


ôn; i 


* [Méc. Anal. 3rd ed. Tome 1, pp. 290-292.] 
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and Ras m(2’ Aorem saie a! a) 
dt òn; Y di Sni  dtdn; (7.) 
-5 m(z vy poy ie) = ôT : 

Sn * oni òn; 


T being here considered as a function of the 6n quantities of the forms y’ and 7, obtained by 
introducing into its definition (4.) the values 
, Ox òx ba 
+o <— +--+ +5n a> Ke. 8. 
=115— 51 Ne Sna 13 San ( ) 


A different proof of this important transformation (3.) is given in the Mécanique Analytique. 


3. The function 7’, being homogeneous of the second dimension with respect to the quan- 


tities 7’, must satisfy the condition 


2T = 2.7 a (9.) 


and since the variation of the same function T' may evidently be expressed as follows, 


stad (oe by +o 5), (10.) 
òn irr 
we see that this variation may be expressed in this other way, 
òT òT 
ST z(a ld S in). (11.) 
If then we put, for abridgement, 
òT òT 
oe? i, (12.) 
and consider 7 (as we may) as a function of the following form, 
T= F (0, Wg, +++ Dgn 71> 22> ++ Nan) (13.) 
we see that 
SF 3; ae 
Su, 72? eee Sty, a (14.) 
and 
sF $87 òF òT 
n a> ' — = -—s 15. 
dn, on Snan ONsn we) 
and therefore that the general equation (3.) may receive this new transformation, 
dw; 5(U-F) 
—= F 16. 
If then we introduce, for abridgement, the following expression H, 
H=F- U= F (w, W2, +e. Wan. M1> 72> +> Nan) — U (M1325 +++ Nan)» (17.) 
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we are conducted to this new manner of presenting the differential equations of motion of a 
system of n points, attracting or repelling one another:* 


dn, 3H. de, 8H. 


ee Oe ea 
dg ae ae ee 
Ge Gta Mk tee (A.) 
Angn 98H | dog, _ _ H 
Ot? TET i By) 


In this view, the problem of mathematical dynamics, fora system of n points, is to integrate a 
system (A.) of 6n ordinary differential equations of the first order, between the 6n variables 
ni, w; and the time t; and the solution of the problem must consist in assigning these 6n variables 
as functions of the time, and of their own initial values, which we may call e;, p;. And all these 
6n functions, or 6n relations to determine them, may be expressed, with perfect generality and 
rigour, by the method of the former Essay, or by the following simplified process. 


Integration of the Equations of Motion, by means of one Principal Function. 


4. If we take the variation of the definite oe 


E =f (2-0; ya = -H) de (18.) 


without varying ¢ or dt, we find, by the Calculus of Variations, 


t 
5s-{ 5S’ . dt, (19.) 
0 
in which 
ate 8H 
S'=2.0—-—H, (20.) 
and thereforet 
òH ôH 
59’ = =3( azo- 5, n). 21. 
ig | (21.) 
that is, by the equations of motion (A.), 
d dw d l 
88’ =x (wo at ae èn) =g . ody; (22.) 
the variation of the integral S is therefore 
5S = È (wdy — pde), (23.) 


* [The canonical form (A.) had already been introduced by Lagrange in the particular case of the equations for 
the variations of the elements in his perturbation theory. Cf. Mé. Anal. 3rd ed. Tome 1, p. 310; also Whittaker, 


Analytical Dynamics (1927), p. 264. n 


t [s= z (E bot wò 5) - -2E io -z èr | 
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(p and e being still initial values,) and it decomposes itself into the following 6n expressions, 
when S is considered as a function of the 6n quantities y;, e;, (involving also the time,) 


bi RBs hig < 88 
w= òn ? Prr Se,’ 
5S 5S 
gag cael 7%, (B.) 
5S 58 
Man Son’ P3n i 


which are evidently forms for the sought integrals of the 6n differential equations of motion 
(A.), containing only one unknown function S. The difficulty of mathematical dynamics is 
therefore reduced to the search and study of this one function S, which may for that reason 
be called the Principat Function of motion of a system. 


This function S was introduced in the first Essay under the form* 
t 
jz Í (T+U)dt, 
0 


the symbols T and U having in this form their recent meanings; and it is worth observing, that 
when S is expressed by this definite integral, the conditions for its variation vanishing (if the 
final and initial coordinates and the time be given) are precisely the differential equations of 
motion (3.), under the forms assigned by Lagrange. The variation of this definite integral S has 
therefore the double property, of giving the differential equations of motion for any trans- 
formed coordinates when the extreme positions are regarded as fixed, and of giving the integrals 
of those differential equations when the extreme positions are treated as varying.} 


5. Although the function S seems to deserve the name here given it of Principal Function, 
as serving to express, in what appears the simplest way, the integrals of the equations of motion, 
and the differential equations themselves; yet the same analysis conducts to other functions, 
which also may be used to express the integrals of the same equations. Thus, if we put 


Q= [,(- x. n'y, +H) de, (24.) 
and take the variation of this integral Q without varying t or dt, we find, by a similar process, 
8Q =X (nda-—edp); (25.) 
so that if we consider Q as a function of the 6n quantities w;, p; and of the time, we shall have 6n 
expressions 
ge ha SARRO: 
N= Em,’ SH Spy? (26.) 


which are other forms for the integrals of the equations of motion (A.), involving the function Q 
instead of S. We might also employ the = 


y= [3.05 w= a i= zf od, (27.) 


* [P. 160.] up [ež 73] - Laff gz or T _ôU D) ay Ja ] 
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which was called the Characteristic Function in the former Essay, and of which, when considered 
as a function of the 6n + 1 quantities n;, e;, H, the variation is ! 

òV =X (wdyn— pde)+t5H. (28.) 
And all these functions S, Q, V, are connected in such a way, that the forms and properties of 
any one may be deduced from those of any other.* 


Investigation of a Pair of Partial Differential Equations of the first Order, 
which the Principal Function must satisfy. 


6. In forming the variation (23.), or the partial differential coefficients (B.), of the Principal 
Function S, the variation of the time was omitted; but it is easy to calculate the coefficient 


“4 corresponding to this variation, since the evident equation | 


St 
dS 88 -èS d 
gives, by (20.), and by (A.), (B.), 
ae yy òH 
ga -2.0 vs =—H. (30.) 


It is evident also that this coefficient, or the quantity — H, is constant, so as not to alter 
during the motion of the system; because the differential equations of motion (A.) give 
dH SH dy H do 
EIRA EA N E 
If, therefore, we attend to the equation (17.), and observe that the function F is necessarily 
rational and integer and homogeneous of the second dimension with respect to the quantities 
w;, we shall perceive that the principal function S must satisfy the two following equations 
between its partial differential coefficients of the first order, which offer the chief means of 
discovering its form: 
5S òS ÒS 5S 
SE (5 i o Sei Ma Mon) = U (715725 +++ Man), 
58 8S 58 58 
Et fe he sey Gigs os 
* (S=V—Ht, Q=—S+X(nw—ep). The process is analogous to that on pp. 174, 175, Mathematical Papers, 
Vol. 1, in which the functions W and T are introduced. See also Routh, Advanced Rigid Dynamics (1905), Art. 487.] 


(C.) 
eon} = U (€,, €g, «2+ gn). 


t [It is now more usual to write S= if (S055 -H) dt, and we have, analogous to (29) and (30), 
ty 
dS 28. os de’ os 0H, es 
digs ch, tE 36 dig ay, (2P Gt Ho) +2P Go Ho 


The two partial differential equations are therefore 
as (0S oes os 
+H SS aaa on 11) =0, 
at Nang” Da aa MY e 
as oS as as 
#(-5 ears tae MEN ERE E PA to) =0, 
even when H is not constant. It is easy to see that these are equivalent to (C.) when H is constant. Jacobi 
(Crelle, xxvul, pp. 97-162) showed that equations (C.) still hold when U contains ¢ explicitly, and the corre- 
sponding equations for the function V are also given there.] 
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Reciprocally, if the form of S be known, the forms of these equations (C.) can be deduced 
from it, by elimination of the quantities e or 7 between the expressions of its partial differential 
coefficients; and thus we can return from the principal function S to the functions F and U, 
and consequently to the expression H, and the equations of motion (A.). 


Analogous remarks apply to the functions Q and V, which must satisfy the partial differential 
equations, 


Sa et ea 


+ (32.) 
Q Q èQ 8Q 82 82 _ 8@ | 
; et (Pas Pas -Pan Sp,’ õp, - 5.) | dp,’ dp,’ =al 
an 
8V 5 ô ) 
P (Sager go tas Mas Man) =H +U (sna aes we 


èV èy èy ) 
P(E Bac’ Gi, Ons ved ean) =H + U (€,, €g; --- ga). 


General Method of improving an approximate Expression for the Principal 
Function in any Problem of Dynamics. 
7. If we separate the principal function S into any two parts, 
S8,+8,=S, (34.) 
and substitute their sum for S in the first equation (C.), the function F, from its rational and 
integer and homogeneous form and dimension, may o an in this new way, 


SE A N =F (Sts g Ja + Mh a) 


+P (mtt (a) ie? | 
ari mmj rin] 
Ea ts (32) 
r(e) -rC 0 
OE 
and since, by (A.) and (B.), 
r (=r o= - ty (38.) 
* [By r (5) is meant the partial derivative of F with respect to a | 
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we easily transform the first equation (C.) to the following, 


dS, a 58, ôS, (eer 58, 
a 7 +U Ms» Nan) — F (ge “Sma tan) +F 8,’ ” Bian =? M19 *** Van] > (D.) 


which gives rigorously 


s= |f- P14 U (4s -N3n) — P (ga Fett on} 


+f F(F... pote tan) dt, 


supposing only that the two parts S,, S,, like the whole principal function S, are chosen so as 
to vanish with the time. 

This general and rigorous transformation offers a general method of improving an approxi- 
mate expression for the principal function S, in any problem of dynamics. For if the part S, be 
such an approximate expression, then the remaining part S, will be small; and the homo- 
geneous function F involving the squares and products of the coefficients of this small part, 
in the second definite integral (E.), will be in general also small, and of a higher order of small- 
ness; we may therefore in general neglect this second definite integral, in passing to a second 
approximation, and may in general improve a first approximate expression S, by adding to 
it the following correction, 

58, 


A8y= | È- SEHU Chs oon) (Sse got mis tan) ba F.) 


in calculating which definite integral we may employ the following approximate forms for the 
integrals of the equations of motion, 
n=- h n=- gh eee Pau ~ es (39.) 
expressing first, by these, the variables n; as functions of the time and of the 6n constants 
€; Pi, and then eliminating, after the integration, the 3n quantities p;, by the same approximate 
forms. And when an improved expression, or second approximate value S,+A8,, for the 
principal function S, has been thus obtained, it may be substituted in like manner for the first 
approximate value S,, so as to obtain a still closer approximation, and the process may be 
repeated indefinitely. 

An analogous process applies to the indefinite improvement of a first approximate expression 
for the function Q or V. 


(E.) 


Rigorous Theory of Perturbations, founded on the Properties of the Disturbing 
Part of the whole Principal Function. 


8. If we separate the expression H (17.) into any two parts of the same kind, 


H,+H,=H, (40.) 
in which 
a H= Fy (01, Dg, +++ Ogn N1 M29 -++ Nan) — U1 (M15 Na» +++ Nan)» (41.) 
an 
H,= F; (0, 02, +++ Man > N1: N23 -++ Nan) — Us (1, N2 +++ Nan)» (42.) 
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the functions F,, Fa, U,, U, being such that 
F,+F,=F, U,+U0,=0U; (43.) 
the differential equations of motion (A.) will take this form, 
dy, 5H, 5H, do; ôH, ôH, 
dt 80, 80,’ dt 8, 8,’ 
and if the part H, and its coefficients be small, they will not differ much from these other 
differential equations, 


dn; 5H, dw; èH, 
d õm dt By,’ 
so that the rigorous integrals of the latter system will be approximate integrals of the former. 
Whenever then, by a proper choice of the predominant term H,, a system of 6n equations such 
as (H.) has been formed and rigorously integrated, giving expressions for the 6n variables 7;, w; 


as functions of the time ¢, and of their own initial values e;, p;, which may be thus denoted: 
i= Pi (t, C15 as +++ Cans P1: Pas +++ Pan)» (44.) 


(H.) 


and 

W; = Pi; (t, 1, Cg, +++ Can» P1 Pa» +++ Pan); (45.) 
the simpler motion thus defined by the rigorous integrals of (H.) may be called the undisturbed 
motion of the proposed system of n points, and the more complex motion expressed by the 
rigorous integrals of (G.) may be called by contrast the disturbed motion of that system; and to 
pass from the one to the other, may be called a Problem of Perturbation. 


9. To accomplish this passage, let us observe that the differential equations of undisturbed 
motion (H.), being of the same form as the original equations (A.), may have their integrals 
similarly expressed, that is, as follows: 


OSA OSA 
w; = ; SS mt aa a I . 
i Sn; Pi 5 ei ( ) 
S, being here the principal function of raame motion, or the definite integral 


considered as a function of the time vale of the vc n> €i- In like manner if we represent 
by S, +S, the whole principal function of disturbed motion, the rigorous integrals of (G.) may 
be expressed by (B.), as eid 

Ra 582 8S, S58, 

-a K. 

“Bag Sy’ PO Beg Be kas 
Comparing the forms (44.) with the second set of equations (I.) for the integrals of undisturbed 
motion, we find that the following relations between the functions ¢;, Sı must be rigorously and 
identically true: ; i 39. 

g 
Mm=ġi (. E1, Cg, +++ Cans — eit Set sie a); (47.) 

and TO BA by (K.), that the integrals of disturbed motion may be put under the following 
forms, 


òS, 55s 85, ) a.) 


n= te (tersen -+ ean Pat Se 2 Pat 3a, 777 Pant Se 


22-2 
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We may therefore calculate rigorously the disturbed variables n; by the rules of undisturbed 
motion (44.), if without altering the time ż, or the initial values e; of those variables, which 
determine the initial configuration, we alter (in general) the initial velocities and directions, by 
adding to the elements p; the following perturbational terms, 
_ ôS, dS, 58, | 
AM= 5 gs a Y ty AP n= 50. * 
a remarkable result, which includes the whole theory of perturbation. We might deduce from 
it the differential coefficients 7}, or the connected quantities w, which determine the disturbed 
directions and velocities of motion at any time t; but a similar reasoning gives at once the 
general expression, 


(M.) 


58, dS, 5S, se N.) 


m= 5+ Ha (reas ear esn: Pit Se, *, Pot Bey’ Pant so 
implying, that after altering the initial velocities and directions or the elements p; as before, by 
the perturbational terms (M.), we may then employ the rules of undisturbed motion (45.) to 
calculate the velocities and directions at the time ż, or the varying quantities w;, if we finally 


apply to these quantities thus calculated the following new corrections for perturbation: 
88, 5S, dS, 


Am =, Aw,==—, ... Aw,=——. 
q 1 z dn» T Nan 


Approximate expressions deduced from the foregoing rigorous Theory. 


(0.) 


10. The foregoing theory gives indeed rigorous expressions for the perturbations, in passing 
from the simpler motion (H.) or (I.) to the more complex motion (G.) or (K.): but it may seem 
that these expressions are of little use, because they involve an unknown disturbing function 
S,, (namely, the perturbational part of the whole principal function S,) and also unknown or 
disturbed coordinates or marks of position »;. However, it was lately shown that whenever 
a first approximate form for the principal function S, such as here the principal function S, of 
undisturbed motion, has been found, the correction S, can in general be assigned, with an 
indefinitely increasing accuracy; and since the perturbations (M.) and (O.) involve the disturbed 
coordinates y; only as they enter into the coefficients of this small disturbing function S,, it is 
evidently permitted to substitute for these coordinates, at first, their undisturbed values, and 
then to correct the results by substituting more accurate expressions. 


11. The function S, of undisturbed motion must satisfy rigorously two partial differential 
equations of the form (C.), namely, 
òS ôS 
yy ts (= eee m 1? Tan) = U4 (7, +++ Nn); 
èS, , p (851 581 

"i ETA 

and therefore, by (D.), the disturbing function S, must satisfy Hecngnaly the following other 
condition: 


d8, òs 5S 
Bs Me ae (stems 
dt ahs: ie A ôn ÒNan 


(P.) 
ee ean) = U4 (645-0 Can); 


21s ota) + (5 pote an) (Q.) 
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and may, on account of the homogeneity and dimension of F, be approximately expressed as 
follows: 


s MA ue BA, 
Sa -f | U, GAs Nan) — F igs: Bi Siren USERE on) a (R.) 
or thus, by (I.), 
t 
S= | {Uatn ++» Nan) — Fa (T1; ees Dyn. N1 ++ nn) (S.) 
that is, by (42.), 
t 
0 


In this expression, H, is given immediately as a function of the varying quantities n;, w, but 
it may be considered in the same order of approximation as a known function of their initial 
values e,, p; and of the time ¢, obtained by substituting for y;, w; their undisturbed values (44.), 
(45.) as functions of those quantities; its variation may therefore be expressed in either of the 
two following ways: 


SHy= 3 ("5889+ 58 a), (48.) 
ôn 
or 
Ei a AAR, 
bH,= 3 ("52 side 529 p) eð (49.) 


Adopting the latter view, and effecting the integration (T.) with respect to the time, by 
treating the elements e;, p; as constant, we are afterwards to substitute for the quantities p,; their 
undisturbed expressions (39.) or (I.), and then we find for the variation of the iia func- 
tion S, the expression * 


ae -H+ E|- fS rk 528 efi Saa), (50.) 


which enables us to transform the perturbational terms (M.), (O.)into the following approximate 
forms: 


8H y, 5 59S: [*8Hs 
An= = 0 a D eese, E dt, (U.) 
and 
828, (dH 
Awm, =}. K ad. V. 
ata vs hat 


containing only functions and quantities which may be regarded as given, by the theory of 
undisturbed motion. 


12. In the same order of approximation, if the variation of the expression (44.) for an un- 
disturbed coordinate n; be thus denoted, 


Sm = BLE (aers a èp), (51.) 
the perturbation of that coordinate may be expressed as follows: 
8n; 
A <= x . ——A $ W. 
7 Sak (W.) 


* [The object is to express S, in terms of ņ and e, or 5S, in terms of n, Sy, e, de.] 
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that is, by (U.), 


tH ô ‘dH, On tH. 
Aiu nef ERa Fy pc cei: ype oe 
W pı J o Se Spa J o bey ~ 8Psn 0 5€gn 
AE ara 
dp, dei dp, ðe, de, oa SPgn Òe 5€3n/) J o OP1 (52.) 
m APIT AA 


(5 5°78, $ ôni 8S aE: > KA sa) ‘oH, 
Bp egne, Spe degn deg BP gn, BC5n) J 0 SP an 
Besides, the identical equation (47.) gives* 


bn Se S, Sm S, y BFS pri 
Se, 5p, de,5e, Sp, de, de, SP5n Oek 5€3y, : 
the expression (52.) may therefore be thus abridged, 
en _ on ‘8H. y ong Su Hy i 
òp J o SaN e (X.) 
òn; ‘OH, dt +.. in dH, dt, 
Se, J o 8p; * Ben ae 


and shows that instead of the rigorous perturbational terms (M.) we may approximately 
employ the following, 


t 
Ap; = ‘Tench dt, k CX.) 
0 


in order to calculate the disturbed configuration at any time ¢ by the rules of undisturbed 
motion, provided that besides thus altering the initial velocities and directions we alter also the 
initial configuration, by the formula 


Ae, = dt. (Z.) 


dH, 
t 
It would not be difficult to calculate, in like manner, approximate expressions for the disturbed 
directions and velocities at any time t; but it is better to resume, in another way, the rigorous 


problem of perturbation. 


Other Rigorous Theory of Perturbation, founded on the properties of the disturbing 
part of the constant of living force, and giving formulæ for the Variation of 
Elements more analogous to those already known. 


13. Suppose that the theory of undisturbed motion has given the 6n constants e;, p; or any 
combinations of these, x4, Kg, -.. Kens aS functions of the 6n variables y;, w; and of the time ż, 
which may be thus denoted: 


Ki = Xr (t, 71> Nes +++ Nans D1» Wes +++ Dgn), (54.) 
and which give reciprocally expressions for the variables ;, w; in terms of these elements and of 
the time, analogous to (44.) and (45.), and capable of being denoted similarly, 


M= $i (t, Kls K2»... Ken)» Oo, = Yi (t, Kl; Kgs... Ken)3 (55.) 


* (Equation (47.) is differentiated partially with respect to e,, the y’s being kept constant.] 
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then, the total differential coefficient of every such element or function x;, taken with respect 
to the time, (both as it enters explicitly and implicitly into the expressions (54.),) must vanish 
in the undisturbed motion; so that, by the differential equations of such motion (H.), the 
following general relation must be rigorously and identically true: 


dK; (3 8H, ôk; ðH ) 


0= (56.) 


a M NSn ba So 

In passing to disturbed motion, if we retain the equation (54.) as a definition of the quantity 
x;, that quantity will no longer be constant, but it will continue to satisfy the inverse relations 
(55.), and may be called, by analogy, a varying element of the motion; and its total differential 
coefficient, taken with respect to the time, may, by the identical equation (56.), and by the 
differential equations of disturbed motion (G.), be rigorously expressed as follows: 
dx; (x òH, dK; =) 


Oink (Si (AX) 


14. This result (A1.) contains the whole theory of the gradual variation of the elements of 
disturbed motion of a system; but it may receive an advantageous transformation, by the sub- 
stitution of the expressions (55.) for the variables y;, w; as functions of the time and of the 
elements; since it will thus conduct to a system of 6n rigorous and ordinary differential equations 
of the first order between those varying elements and the time. Expressing, therefore, the 
quantity H, as a function of these latter variables, its variation 5H, takes this new form, 
dH, ôH, 
e in Bes 


and gives, by comparison with the form (48.), and by (54.), 
dH, . ôH, x. ôH, ôH, 8x | 


8H, = 23x + ——*8t, (57.) 


Sa. be th, tas Be on eink 
and thus the general equation (A1.) is transformed to the following, 
| dey tn 3H, wah, ôH, ” 
Fae. | fel eR ae st he (B .) 
in which | 
(8K; 8K, Òr; Öre), h 
T (Ge get Seg): (C me 


so that it only remains to eliminate the variables y, w from the expressions of these latter 
coefficients. Now it is remarkable that this elimination removes the symbol ¢ also, and leaves 
the coefficients a; , expressed as functions of the elements «x alone, not explicitly involving the 
time. This general theorem of dynamics, which is, perhaps, a little more extensive than the 
analogous results discovered by Lagrange and by Poisson,* since it does not limit the disturbing 
terms in the differential equations of motion to depend on the configuration only, may be 
investigated in the following way. 


* (Lagrange, “Sur la théorie générale de la variation des constantes arbitraires,” Mémoires de Ul’ Institut (1808), 


p. 288; Poisson, “Sur la variation des constantes arbitraires dans les questions de mécanique,” Journal de l École 
Polyt. Tome vin (1809), pp. 288, 289. See reference to Lovett, p. 160.] . 
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15. The sign of summation X in (C1.), like the same sign in those other analogous equations 
in which it has already occurred without an index in this Essay, refers not to the expressed 
indices, such as here #, s, in the quantity to be summed, but to an index which is not expressed, 
and which may be here called r; so that if we introduce for greater clearness this variable index 
and its limits, the expression (C1.) becomes 

3n (OK; 5K, Örk; Ôk 
aam Zat (5p ba; Be S) P 
and its total differential coefficient, taken with respect to the time, may be separated into the 
two following parts, 
d Ly 3n (z d 8k, òk, d 5) 


di? Ze | Sn, dda, Sy. dl 5a, 


60. 
5 Ae dòk; òk; d a), wis 


1\3o0,dt8n, 5a,dtdn,)’ 
which we shall proceed to calculate separately, and then to add them together. By the definition 
(54.), and the differential equations of disturbed motion (G.), 
d òr SK; 3n | ôk; (3H, èH) ô; (8H, dH, 
diso, tm, W1 i (5 tm, 2) - 50,50, (5 ny alt 
in which, by the identical equation (56.), 
8k; ò a 3n/dx;5H, 5x; 5H, 
Sib, ~~ Bas20 (Spe Sme Be, By) (62) 


(61.) 


dtdm, da, ™1 
we have therefore 


d Bx, yan = (glia 8%, SHa Ski 8Ha dx; SH, Sr; 8H, ‘), ie 
dt8o0, ~“!1\8n, 80,80, 50,50, 5y, Òu ÒN ÒD, 57, 50,50, 
and aie may be found from this, by merely changing 7 to s: so that* 
y 2n (òx; d xy _ By dB 
1\5n,dida, Sn, d S0, 
in axial (2s Sx; Sr; Bk, ) Bez 8%, Bey Sh \3Hy | gy 
“ew 11 | (y, 50,30, Sn, 50,50,) Su * \8n,57,50, 57,57,80,) 5m, | CH 
+ (Bee Bre Bardet) PH (ecb) 
ôn, dm, ÒN, do, ôn, òw, (55 ÒN, on, fates 
and similarly, 
mn (r, d Bq; _ Be: dB, 
1\ 5m, dt Sn, 5w, di m, 
is gifi Bx; Sr; Bx, jz (2 Sky öx, Si \ M| ig 
edi ((50,By,80, 80, Snabn) Bo, * (8 3a, Be,5a,by) Bae | (8) 
(BBs Sse) BH, (Be dau _ tu By) SA 
80,5, ÒT, ÒN, 50,5, (T, S, ÒT, ÒT, 5n,,5n,) ° 
* (In (64.) and (65.) the coefficients of 52. and of gt are easily seen to be zero, Also if we interchange 


2H, 
the suffices in the coefficient of ie Sa, in (64.) the term cancels against the corresponding term of (65.).] 
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Adding, therefore, the two last expressions, and making the reductions which present them- 
selves, we find, by (60.),* 


ou = Sind (4 = + BY 7 4), (D!) 
in which 
P a 3n (5 Bx; SK; OK, + OK Ox, dK, Òk; 
1\$n,50,50, 57,50,50, 50,50,5n, ÒT, ÒT, ÒN, 


(66.) 


pY RA (5: Ox; — 8K; T + Oke 8x, dK, 8x; ); 
ôw, nun, 50, 87.5, Ôn, 89,50, ÒN, 9,80, 


and since this general form (D1.) for Sa s contains no term independent of the disturbing 


quantities Se $ e, it is easy to infer from it the important consequence already mentioned, 


namely, that the coefficients a; ,, in the differentials (B+.) of the elements, may be expressed as 
functions of those elements alone, not explicitly involving the time. 


It is evident also, that these coefficients a, , have the property 


as, iT — lis» (67.) 
and 
Ai, r 0; (68. ) 
the term proportional to ——* | ? disappears therefore from the expression (B1.) for = ; and the term 


DA 

SHy 8H, . SHydx; 

Bic, 8" Bag Oy. aE 
destroys the term 

SH, SH, , Beds, 

8x, &** 8x; dx, dt’ 


hima Seis st 1 
yeaa (EL) 
or 
dH, _èH,, i 
s O87 (E.) 


that is, in taking the first total differential coefficient of the disturbing expression H, with 
respect to the time, the elements may be treated as constant. 


Simplification of the differential equations which determine these gradually varying 
elements, in any problem of Perturbation; and Integration of the simplified 
equations by means of certain Functions of Elements. 


16. The most natural choice of these elements is that which makes them correspond, in 


* [It follows that in the undisturbed motion the a’s are constant, which is Poisson’s theorem (Journal de l’ École 
Polyt. Tome vrn (1809), pp. 281, 282).] 

t [This is most easily seen as follows. If a; ,, when expressed in terms of the xs, contains ¢ explicitly, it would 
not remain constant during the undisturbed motion, which would contradict (D1.).] 


HMPII 23 
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undisturbed motion, to the initial quantities e;, p; . These quantities, by the differential equations 
(H.), may be expressed in undisturbed motion as follows, 


t H ‘6H 


and if we suppose them found, by elimination, under the forms 


= Ni +D; (t,71,72, eee Ngn» @1> Bo, eee ici (70 ) 


Di = 0; + Ey (t, 715 Nas +++ Nan» D1» D2, +++ Oan), 
it is easy to see that the following equations must be rigorously and identically true,* for all 
values of n;, Oi, 

pe ala cH Hon! (71.) 


o=; (0,71, %, eee Ngn> Py; Do; eee wgn): 


When, therefore, in passing to disturbed motion, we establish the equations of definition, 


ki=nNi +Ò; (t, N1, Nes ++ Nans D1, Da, +++ Dzn), | (72.) 


À= @ + Fi (t, 71; N25 +++ Nan» D1» Wg, +++ Dyn); 
introducing 6n varying elements «;, A;, of which the set A; would have been represented in our 
recent notation as follows: 

Ài = Kon43 (73.) 


we see that all the partial differential coefficients of the forms we oe i i é = , vanish when 
dn, òm, dn,’ 80, 
t=0, except the following: 
oi ea (74.) 
òn; 8m; 
and, therefore, that when ¢is made = 0, in the coefficients a; ,, (59.), all those coefficients vanish, 
except the following: . 
Dy, 3n+r = l, Anr, = — l. (75.) 


But it has been proved that these coefficients a;,, when expressed as functions of the 
elements, do not contain the time explicitly; and the supposition ¢=0 introduces no relation 
between those 6n elements «;, A;, which still remain independent: the coefficients &; ,, therefore, 
could not acquire the values 1, 0, — 1, by the supposition ¢ =0, unless they had those values 
constantly, and independently of that supposition. The differential equations of. the forms 
(B1.) may therefore be expressed, for the present system of varying elements, in the following 
simpler way: ; 

dk; H, ddA;  dHy. 

Pa ee war A, 
and an easy verification of these expressions is offered by the formula (E*.), which takes now 
this form, 


(Gt) 


pi dx ôH, 4 nl 


i ee ae bees 1 
Se dit OA di (H?.) 


* [Putting t=0 in (70.) we have ®; (0, e1, ... egn» P1» --- Pan) =Q for all values of the e’s and p’s. Hence (71.) 
follows.] 
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17. The initial values of the varying elements «;, A; are evidently e;, p;, by the definitions 
(72.), and by the identical equations (71.); the problem of integrating rigorously the equations 
of disturbed motion (G.), between the variables n;, w; and the time, or of determining these 
variables as functions of the time and of their own initial values e;, p;, is therefore rigorously 
transformed into the problem of integrating the equations (G".), or of determining the 6n 
elements «;, A; as functions of the time and of the same initial values. The chief advantage of 
this transformation is, that if the perturbations be small, the new variables (namely, the 
elements,) alter but little: and that, since the new differential equations are of the same form 
as the old, they may be integrated by a similar method. Considering, therefore, the definite 
integral 


t 
E= f € A= Hs) dt, (76.) 


as a function of the time and of the 6n quantities x, , Kg, ... Kgn» 15 Cg) +++ gn» and observing that 
its variation, taken with respect to the latter quantities, may be shown by a process similar to 
that of the fourth number of this Essay to be 

SE = X (ôx — põe), (I*;) 
we find that the rigorous integrals of the differential equations (G!.) may be expressed in the 


following manner: 
òE SE 


——-_— p T mip =a 1 
À; SK,’ Pi Se; ? (K .) 


in which there enters only one unknown function of elements E, to the search and study of which 
single function the problem of perturbation is reduced by this new method. 


We might also have put 


t 
ô=] (-2. +H) d, (17.) 
0 òk 


and have considered this definite integral C as a function of the time and of the 6n quantities 
Ài, pi; and then we should have found the following other forms for the integrals of the differential 
equations of varying elements, 


(L*.) 


And each of these functions of elements, C and E, must satisfy a certain partial differential 
equation, analogous to the first equation of each pair mentioned in the sixth number of this 
Essay, and deduced on similar principles.* 


18. Thus, it is evident, by the form of the function E, and by the equations (K".), (G1.), and 
(76.), that the partial differential coefficient of this function, taken with respect to the time, is 


BH _dE_y BE de _ 
| St dt ~'8kdt 

and therefore that if we separate this function Æ into any two parts 
j E + E= E, (N1.) 


— H,; (M".) 


* [0 and E are the generating functions for the contact transformations from p, e to A, «.] 
3 23-2 
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and if, for greater clearness, we put the expression H, under the form 


HF gts Fh (0,003, Begs oan s My s Ay o> Agy)s (O1.) 
we shall have rigorously the partial differential equation 
ôE, ôF, ôE, ôE, ôF, ôE, y 
0 r St 24 Halt rs Kgn> © Sk, Trigg ng tg}: (P cq 


which gives, approximately, by (G.) and (K1.), when the part Æ, is small, and when we neglect 
the squares and products of its partial differential coefficients,* 


dE, SE, SE, 8E, 
0= Get apt Halter Kan: aot na) 


Hence, in the same order of approximation, if the part Æ, like the whole function Æ, be chosen 
so as to vanish with the time, we shall have 


z ôE ôE, ôF ; i 
E,= — -f ee +H, (tr Kon? Be 15) fat (R1.) 


and thus a first approximate expression E, can be successively and indefinitely corrected. 


(Q?.) 


Bega 


Again, by (L1.) and (G+.), and by the definition (77.), 


80 dO 4 80d y : 
ae ee Sy ge Ee (S1.) 


the function C must therefore satisfy rigorously the partial differential equation, 


ôC 50 òC 
Ha (* 5A,” saaga? ay Dan) (T!) 
and if we put 
C=C,+C,, (0) 
and suppose that the part C, is small, then the rigorous equation 
òC 80, 80, 400; 80, +002 j 
t =H,(t, TA ns Wy We don) care 
becomes approximately, by (G1.) and (L?.), 
dC, 80, èC, 8c, . 
net lt + i Bae s Ags . T VOR , (W1.) 
and gives by integration, 
soft gr 8; Cı ôC, r 
o= f d- i t+ H(t, DE EA Aon) bd (X1) 


the parts C, and C, being supposed to vanish separately when t= 0, like the whole function of 
elements C. 


And to obtain such a first anh ‚E, or C,, to either of these two functions of 
elements Æ, C, we may change, in the definitions (76.), (77.), the varying elements x, A, to their 


* [The deduction is the same as that of (D.), p. 170.] 
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initial values e, p, and then eliminate one set of these initial values by the corresponding set of 
the following approximate equations, deduced from the formule (G1.): 


‘SH, 
se oTe dt; yi, 
Ki e+ ‘ Sp, dt; ( ) 
and 
a ‘SH, . 
=p- f Be dt. (Z}.) 


It is easy also to see that these two functions of elements C and E are connected with each 
other,* and with the disturbing function S,, so that the form of any one may be deduced from 
that of any other, when the function S, of undisturbed motion is known. 


Analogous formule for the motion of a Single Point. 


19. Our general method in dynamics, though intended chiefly for the study of attracting 
and repelling systems, is not confined to such, but may be used in all questions to which the law 
of living forces applies. And all the analysis of this Essay, but especially the theory of per- 
turbations, may usefully be illustrated by the following analogous reasonings and results 
respecting the motion of a single point. 

Imagine then such a point, having for its three rectangular coordinates x, y, z, and moving in 
an orbit determined by three ordinary differential equations of the second order of forms 
analogous to the equations (2.), namely, 

” ôU d ” rae dU Fe ” ô i 


t= Sz > j dy ee ~ 5 
U being any given function of the coordinates not expressly involving the time: and let us 
establish the following definition, analogous to (4.), 

T =} (x+y? +z"), . (79.) 
x’, y’, 2’ being the first, and x”, y”,z” being the second differential coefficients of the coordinates, 
considered as functions of the time ¢. If we express, for greater generality or facility, the 
rectangular coordinates x, y, z as functions of three other marks of position 7, , 42, 43, T will 
become a homogeneous function of the second dimension of their first differential coefficients 
"i> 12, Nng taken with respect to the time; and if we put, for abridgement, 
òT òT òT 


(78.) 


a ont ws tae T oe On’ (80.) 
T may be considered also as a function of the form 
T = F (1, 02, 3,715 25%) (81.) 


which will be homogeneous of the second dimension with respect to ©, Dz, Dz. We may also put, 
for abridgement, 
F (@,, D2, 03, M1: N2: N3) — U (m5 92593) =H; (82.) 


and then, instead of the three differential equations of the second order (78.), we may employ 
* [C=E (nA—ep)—#.] 
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the six following of the first order, analogous to the equations (A.), and obtained by a similar 
reasoning, 
dm 3H dm_ ôH dm_ 8H 
dt ôm, dt ôm, dt 50,” 
do, ôH da, ôH dos SH 
dt ~~ 8m,’ dt ~~ O S 
20. The rigorous integrals of these six differential equations may be expressed under the 
following forms, analogous to (B.), 


(83.) 


ye ARo ode goe AN 
1 Sn,’ 2 ine’ a ons’ 
1 2 3 
5S 5S 5S ot) 


P= ~ Be,’ P= i! Ps= ruda? 
in which e; , €s, €3, P1» Pa» Pa are the initial values, or values at the time 0, of 11, Na» 73, Drs Da, D3; 
and S is the definite integral 
‘/ dH 8H 8H 

s=| (055 + 55, t 55,2) at (85.) 
considered as a function of 7, , 72,3, €&1» 2, €g and t. The quantity H does not change in the course 
of the motion, and the function S must satisfy the following pair of partial differential equations 
of the first order, analogous to the pair (C.), 
òs ôS 8S ös 
aed (ae eee ; ; )=u > ; ; 
St Sm’ dq’ dns M1> N2» N3 (71; Na» Ns) 


ôS 5S 8S 88 

wT (= Se,’ Se,’ Cis êz, es) =U le, ez, €z). 
This important function S, which may be called the principal function of the motion, may 
hence be rigorously expressed under the following form, obtained by reasonings analogous to 
those of the seventh number of this Essay: 


(86.) 


t( $88 8S, 5S, ôS 

S=8. +f [- +U , , -F(S, <, <, , ’ )jar 
at) 1a FU h Me» a) Sm’ Sha’ Smp? 1? 3 
t (8S 88, OR. 88, 88 88 
P(E E a <—- =z, , , ) äs; 
J, õnı Ön Sna Sna’ Sna Snp’ "H "2 9 

S, being any arbitrary function of the same quantities 7, , n2, 73, €, €z, €z, é, 80 chosen as to vanish 
with the time. And if this arbitrary function S, be chosen so as to be a first approximate value 
of the principal function S, we may neglect, in a second approximation, the second definite 
integral in (87.). 

21. A first approximation of this kind can be obtained, whenever, by separating the 
expression H, (82.), into a predominant and a smaller part, H, and H,, and by neglecting the 
part H,, we have changed the differential equations (83.) to others, namely, 

dn 8H, ty 8H, — dm _ 8, 
dt 80,’ dt 8m,’ dt 8a,’ 
do,__ 5H, do, 5H, do, 8H, 

Sna 


(87.) 


(88.) 


di 3m? de Big’ a 
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and have succeeded in integrating rigorously these simplified equations, belonging to a simpler 
motion, which may be called the undisturbed motion of the point. For the principal function of 
such undisturbed motion, namely, the definite integral 


*/ S68, òH, 6H, 
S,= (aa. +O ELA ttmg) dt (89.) 


considered as a function of 71, N2, 73, 1 €z, €z, t, Will then be an approximate value for the original 
function of disturbed motion S, which original function corresponds to the more complex 
differential equations, 


dn,_8H, 3H, — dyg_8H, SH, —dyg_8Hy ôH, 


dt 80, 80,’ dt 8m, Sm,’ dt Sw, dm,’ 
do, èH, SH, dm, __ SH, ôH, do, 5H, èH, 


dt Sm m dt Bg Bg’ dt Og 
The function S, of undisturbed motion must satisfy a pair of partial differential equations of 
the first order, analogous to the pair (86.); and the integrals of undisturbed motion may be 
represented thus, 


(90.) 


CE oe OE ENE. REE GN 
1 87,’ 2 Sne’ 3 873” pa N 
88, 8S, 89, l 
aTe My P2= ye. pimi 


while the integrals of disturbed motion may be expressed with equal rigour under the following 
analogous forms, 


È D gp p8 | BSy 
t èm òm’ * 82 89’ > Sns Öna’ | (92.) 
- ÒS, _ ÒS S, 88, 88, 89, i 
a ee i oa e E Se, Seg 


if S, denote the rigorous correction of S,, or the disturbing part of the whole principal function 
S. And by the foregoing general theory of approximation, this disturbing part or function S, 
may be approximately represented by the definite integral (T.), 


t 
| 0 
in calculating which definite integral the equations (91.) may be employed. 


22. If the integrals of undisturbed motion (91.) have given 


11 = 1 (t, €1, €z, 3, Pr» Po» P3), 
X Na = $2 (t, e1, ĉ2, €33 Pi» Pa» P3), (94.) 


na = 3 (t, €1, 9, &3, P15 P2» P3)» 
and 


Dy = Wy (t, 6&1, e2, &3, P1, Pa» P3)» 
Ty = Yo (t, 1, Cg, 3, P1» Pa» P3), (95.) 
Dg = ps (t, ĉis 9, &3, Pi» Po» Ps)» 
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then the integrals of disturbed motion (92.) may be rigorously transformed as follows, 


ie 58, dS, d8, 
m=d1(t €1, 2, &3, Pits Pot se Pta) 


= 582 òS 5S, 
m= (he Cf» Pits Pats ptg) 
n= $s (i €1, 2, 3, Puhipi Paty? Ps tSo 


(96.) 


and 
5S 
m5, +44 (i ĉi» 2, Cg, Pıt 


LERLE) 
se’ P2 Se,” P3 Se, )” 


553, 88 5S, 581) fon 


i sei ah a ——— 
m= 5+ ba (' eis £2, &3, Pit so Pat so Pa + Ses 


pa dS. +a (i ea. € re òS, + 5S, ‘he 4 
@3= One ON” ĉis 2, C3, Pi Se,’ Pe Se,’ Ps Se, )” 
S, being here the rigorous disturbing function. And the perturbations of position, at any time t, 


may be approximately expressed by the following formula, 


hn, 30 [3H gy 9m ("Bll qp èm [Eag 
1 8e,J 0 8p; Òe J 0 Spe Sez J 0 Ops (08, 
_ Sm fH, gy 5m S'S, y Sm f'E | 
8p, J o Se Spe Jo Ses Sp3Jo Se, 


together with two similar formule for the perturbations of the two other coordinates, or marks 
of position 72,73. In these formule, the coordinates and H, are supposed to be expressed, by the 
theory of undisturbed motion, as functions of the time t, and of the constants e4 , €z, €3, 1, P2» Ps- 
23. Again, if the integrals of undisturbed motion have given, by elimination, expressions 
for these constants, of the forms 
€y= 1+, (t, 21, N2» N3 D1, D2, Ws); 
€2= Na+ Do (t, 11, N2; Ns» M1, Me, =| (99.) 


€3= 3+; (t, i> Ne> 13> M1, We; Wz), 
and 
Pi=9,+ Filt, 11, 2, 13, Dis Da, | 


Pa=0 + V's (t, 215 Ne, Ng G1» Be, W3), 

P3= 03+ Vs (t, 11, Nas 13, M1, Dos ay: 
and if, for disturbed motion, we establish the definitions 

Ky=1 +, (t, m, N2 Ng Tis Da, = 


(100.) 


Ko=Nna +P (t, 11, N2, 13, D1, De, W3), (101.) 


K3 =N + Ds (t, N1> Ne, Ng, Wis Da, Ws), 
and 
A,=0,+ F(t, 71, M2, 03, Dis De; @s); | 
L 


Ag= 0 + Fa (t, 1, Ne; 13, Vis We, Ws); (102.) 
Às = Wg + Ys (t, 11, Ne, 13, D1, De, W3);! 
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we shall have, for such disturbed motion, the following rigorous equations, of the forms (94.) 
and (95.), 
1 = 4 (t, Kis Ka, Kg, Az, Àg, Ag); 


Ne= $a (t, Ky, Ka, Kgs Ay, Ag, As); (103.) 
13 = $3 (t, Ky, Ka, Kg, Ay, Ags As); 

and T,=} (t, Ky, K2, Kg, Az, Ag, Ag); 
y= Yo (t, Ky, Ko, Kg, Az, Ag, Ag), (104.) 


Ws = Ys (t, Ky, Ko, Kg, Ay, Ag, Ag); 
and may call the quantities x, , xa, «3, Ay, Ag, Ag the 6 varying elements of the motion. To determine 
these six varying elements, we may employ the six following rigorous equations in ordinary 
differentials of the first order, in which H, is supposed to have been expressed by (103.) and 
(104.) as a function of the elements and of the time: 

dx, H, dx, 5H, dx; òH, 
BD,’ de Pa T 
dà _ _ dH, da, dH, dr, dH, 


yy 8K,” da Sa a SK,” 


and the rigorous integrals of these 6 equations may be expressed in the following manner, 


(105.) 


oF = LH _oH# 
1 8K? ey eT Big’ 
SH SE SH (106.) 


P= me P2= $e,’ ss a 


the constants ¢,, €z, €z, P1» P2, Pz retaining their recent meanings, and being therefore the initial 
values of the elements x, , Ka, K3, Ài, Àg, Àg ; while the function Æ, which may be called the function 
of elements, because its form determines the ae of their variations, is the definite integral 

E= af: (at rz a Ha! -H 2), (107.) 
considered as depending on «4, Kg, K3, €,, ĉo, €g and t. The integrals of the equations (105.) may 
also be expressed in this other way, 


HO Lda OO aaa = 
i ees i 7 eS 8 s= “Tay. à 
gig Y 5A, 5A, EIN 
pulpal O Ss 
<p tye Pi oe epea TNR oe 
ing the de integr. 
C=- hdd seal i a Ja (109.) 
a. a am iA. 78) + 


regarded as a function of À; , Àz, mi , Pi» P2» Pz and t: and it is easy to prove that each of these two 
functions of elements, C and E, must satisfy a partial differential equation of the first order,* 
* [The partial differential equations are 


ac ec oH oF 
At tals an’ r)=0, Th Hy (t, Gx )= 0.) 


HMPII l i 24 


www.rcin.org.pl 


186 Ill. GENERAL METHOD IN DYNAMICS [23, 24 


which can be previously assigned, and which may assist in discovering the forms of these two 
functions, and especially in improving an approximate expression for either. All these results 
for the motion of a single point are analogous to the results already deduced in this Essay, for 


an attracting or repelling system. 


Mathematical Example, suggested by the motion of Projectiles. 


24. Ifthe three marks of position 7, , 2,73 0f the moving point are the rectangular coordinates 
themselves, and if the function U has the form 


| U = — ghs — $ {u? (ni +12) + m3}, (110.) 
g, p, v being constants; then the expression 
H=} (aj + w3 + w3) + gns + $ {u (mi +13) + ng} (111.) 


is that which must be substituted in the general forms (83.), in order to form the 6 differential 
equations of motion of the first order, namely, 


—=— = ea = —=—_ = 
at's ” iy Ring Hina 
dw dw dw f 
u 7 Te ae — pn; a —g—v'ns. 
These differential equations have for their rigorous integrals the six following, 
=e cos yt +7 sin pt, 
Ne = es cos pt +7 sin ut (113.) 
Nz = £z COS vt +P3sin vt— vers vt, 
and W = Pp; Cos ut — we, sin pt, 
Wa =P. COS pt — ue Sin pt, 
2 = P2 COS pt — peg SIN pt (114.) 


Ws = Pz COS vt—(ve,+2) sin ot 


e1, €; €3, Py, Po, Pa being still the initial values of 71, n2, Nng, W1, Da, Wg. 


Employing these rigorous integral equations to calculate the function S, that is, by (85.) 
and (110.), (111.), the definite integral 


t 2 2 2 
s= (ee u) dt, (115.) 
we find L(af-+ 8+ m8) =1 | pt+ p+ pitate (veat) | 
+E {pit pi — py? (ef + €3)} cos 2ut — du (e,p, + ep.) sin 2ut (116.) 
g\* Ha g 
EAE e r | 
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2 
and v= -4 {ptr pbe i+ itere) | | 
+4{pi+pi— p’ (ei +e3)} cos 2ut— dy (ey + ezpa) sin Zut (117.) 
g\° Nas 
+4 pt- (vea + 2) | cos 2w- 4 (vea +2) pasin 2 | 


and therefore, 
S=% +{pi+pi- (+e) Te 


+{pt- (rat! o) E 8) em | 


In order, however, to express this function S, as supposed by our general method, in terms 
of the final and initial coordinates and of the time, we must employ the analogous expressions 
for the constants p,, Pa, Pa, deduced from the integrals (113.), namely, the following: 


p _ Hy — pe Cos pt 
A sinpt =’ 
p __ Hla — pecos pt 
” sinpt =’ (119.) 
apt g 
ist) (ves +8) cos ve 


atoe T a 


(118.) 


and then we find 
g= gt pe. (m — 61)? + (12 — 22)? Fod (13 — e3)? 
22 2° tan pt 2° tanvt 
(120.) 


-pme tme) tan! —» (no + 2) (a+ 2.) tan 5, 


This principal function S satisfies the following pair of partial differential equations of the 
first order, of the kind (86.), 


eta (ge) +(5) HEY |= = Ms -E tend) — n, 


ôS 1(/88 ôs 5S\?) | We; % 
a +a (a) +a) + (66) |= - — Ft) - 
and if its form had been previously found, by the help of this pair, or in any other way, the 


integrals of the equations of motion might (by our general method) have been deduced from it, under 
the forms, 


(121.) 


a= BS u(n- e)ootanpt pe, tant, 


5S 
= 5, =H (ney) cotan pt— pestan fa, (122.) 


I) tan” 
tan 5 


v 


èS 
= 5p, 7” (1a ea) cotan vt- (ve, + 


24-2 
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òS t 
and e - 5g, =H (me) cotan pt + wn, tan, 
ôs 
P= ~ 5, TH (nae) cotan pt + patan fy, (123.) 


58 t 
P= -5e T” (me) cotan vt+ (rng +2) tan: 


the last of these two sets of equations coinciding with the set (119.), or (113.), and conducting, 
when combined with the first set, (122.), to the other former set of integrals, (114.). 


25. Suppose now, to illustrate the theory of perturbation, that the constants p, v are small, 
and that, after separating the expression (111.) for H into the two parts, 


H, =} (oj + 03+ 05) +903, (124.) 


A=} {u° (ni + 3) + vni} (125.) 


we suppress at first the small part H,, and so form, by (88.), these other and simpler differential 
equations of a motion which we shall call wndisturbed: 


dm_, Mo, % 


and 


(126.) 


These new equations have for their rigorous integrals, of the forms (94.) and (95.), 


$ Mm=Cyt+Pyt, Na=ea+ Pat, Ng = lz + Pzt- hg, (127.) 
an 


Tı=Pı; G=P2, Wz =Pz— Jt; (128.) 
and the principal function S, of the same undisturbed motion is, by (89.), 


t 2 2 2 
we (H gn) at 


t 2 2 2 
i i Í (2 +0 H ot ge) dt (129.) 
0 f 


2 2 2 
= (Atte E ges) ; att IPs {2 $ 49743, 
or finally, by (127.), 


ki (1—6)? + (ne rs e2)? + (n3— e3)? _ hot (no +e) -hg P. (130.) 


This function satisfies, as it ought, the following pair of partial differential equations, 
3a TO E E fara 
ôt 2 òn dn» ôns | (131.) 


aS, HS (8) Sag 
ôt 2 |\ 8e JA PAS B adii 
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And if by the help of this pair, or in any other way, the form (130.) of this principal function 8, 
had been found, the integral equations (127.) and (128.) might have been deduced from it, by 


our general method, as follows: 


5S. — 

e ds 

pomt, (132.) 
mmga pgh, 

and 

p=- pnrt, 
fait -j= noe, (133.) 
p= -FBT Bot 


the latter of these two sets coinciding with (127.), and the former set conducting to (128.). 


26. Returning now from this simpler motion to the more complex motion first mentioned, 
and denoting by S, the disturbing part or function which must be added to S, in order to make 
up the whole principal function S of that more complex motion; we have, by applying our 
general method, the following rigorous expression for this disturbing function, 


88)? (382)? , (38 
Hatt] 3 (in) * in) * Bre) 12 
=- Í, on Ons dns | 
in which we may, approximately, neglect the second definite integral, and calculate the first 
by the help of the equations of undisturbed motion. In this manner we find, approximately, 


by (125.), (127.), 


(134.) 


2 2 
-H,= — 5 {(e1 + p14)? + (ez + Pt)? <3 (e3 + pt — $gt*)?, (135.) 
and therefore, by integration,* 
Sa = — 4 {u? (e2 +e) + v3} t— 4 {u? (e py + eapo) + vespa} t? | (136.) 
—}{u? (pi + p3) + v? (p3 —gea)} B+ Sop — dorg P, 


* [The approximate S, of (136.) is primarily a solution of the partial differential equation 
8Sa , 38; 554 | 58; 55s, 55; 355, 7 =O; 
8 8m na Ne dns Öns 
It must be a complete integral, the arbitrary constants being ¢,, £2, ¢, and an additive constant. The method 
employed by Hamilton in such cages is as follows. The value of S, is found by integrating along the path 


S, is then a function of t, e1, €2, €g» P1» P2» Pae Replacing p, by — i , etc., S now satisfies the partial differential 


equation because the part arising from the differentiation of the p’s vanishes on account of relations such as 


ap, 55; be , 38: Bp 4 88 BP _o, 
Tei in "Oh 5n2 òns ns ; 
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3 
or, by (133.), S= -p (ni + erm ++i t eant e) | 
j (137.) 
v 
A + erpa ted do (mete) + age): 


the error being of the fourth order, with respect to the small quantities u, v. And neglecting this 
small error, we can deduce, by our general method, approximate forms for the integrals of the 
equations of disturbed motion, from the corrected function S, + S,, as follows: 


_ 88, 88, M-e _ pt 
B= in cd a hes n-a (ns +46), 
ôS: Sa Nz ês 
= sura m + łe), 138. 
2 Sna Sna i salt + 2 ( ) 
_ 88; 5S, _ p 
aeei =n (n+ $e + 
~ 8g Sna — gt A s 3 490 ii 
and 
e AE A Mm 41 pet 
Pi aea < z (4+ 3), 
_ 38, 9 n-e a 
P2= des iter dat | + (ea + $2), (139.) 
ôS, ôS 
p= — So 952 S gtt (eat inat dat) 
3 3 


or, in the same order of approximation, 
y=, +pyt— fut (e1 + $748), 
Na = lg + Pot — Zyl (ea + Pal), (140.) 
Na = C3 + Pat — gt — pvt? (e3 + $pat — vag"), 
and 
w, = Py — pt (e, + 3941), 
Wy = Po— Ht (ez + Pat), (141.) 
wz = Pz — Jt — v*t (e3 + 4 P3t— bgi). 
Accordingly, if we develope the rigorous integrals of disturbed motion, (113.) and (114.), as far 
as the squares (inclusive) of the small quantities u and v, we are conducted to these approximate 


integrals; and if we develope the rigorous expression (120.) for the principal function of such 


motion, to the same degree of accuracy, we obtain the sum of the two expressions (130.) and 
(137.). 


27. To illustrate still further, in the present example, our general method of successive 
approximation, let S, denote the small unknown correction of the approximate expression 
(137.), so that we shall now have, rigorously, for the present disturbed motion, 


S=8,+8,+S8s, (142.) 
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S, and S, being here determined rigorously by (130.) and (137.). Then, substituting S, + S, for 
S, in the general transformation (87.), we find, rigorously, in the present question, * 


#1(/d8,\? [98 dS, 
a a * (in) 1% 
” 0 2 (\dn, On òns 
ee >) (=) (= y | 
+} 54x) +(K2) + dt: 
f 02 (la SNe dns 
and if we neglect only terms of the eighth and higher dimensions with respect to the small 
quantities u, v, we may confine ourselves to the first of these two definite integrals, and may 


employ, in calculating it, the approximate expressions (140.) for the coordinates of disturbed 
motion. In this manner we obtain the very approximate expression, 


at 
S= — = Ka {(m1 + $61)? + (nq + $e2)°} dt 


(143.) 


f e mtteatigea 


=e E Ant + Tye, +408 + Anh + Tey + Aet ) 


p 
Viva 


vigt (144.) 
240 (18+ &s)— -10320 


~~ 360 vit" (ant P 73 e3 + 4e2) — 
-E t+ Hn Hime + ef +: 3 + Hne +e?) 
+ ee i 316g 
-gag (18 +1828 +63) — ‘a Ns +s) — 795760: 


which is accordingly the sum of the terms of the fourth and sixth dimensions in the development 
of the rigorous expression (120.), and gives, by our general method, correspondingly approxi- 
mate expressions for the integrals of disturbed motion, under the forms 


_ 88, , 88, 3S, 
+..-* 


w 2p 
vi bn õn Òn 
_ oS; 88, 58, 
243 (145.) 
~ 8g 82 Sno” 
i _ ôS 58, 58, 
rg ôns ôns ôns Na’ 
* [S, and S, are rigorously determined by 
71445 a) =-H, (i) 
and 4 a +554 ma =- H}. (ii) 
1 
In (143) 2 ¢ means qth i where S=8, +S, +S, and so becomes 
O85 , y (29, 3Sa , 395) AS aS, BY. n 
a +2 (3 ata) Fete (a j- “EG a)! 0 (iii) 
We see at once that the sum of three equations (i), (ii) and (iii) gives © e+e (E) = —H,—H,.] 
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ye 8S, dS, ôS; 
and Pno Bark Ta? 
_ 38; 8S 3 
P= — Se, Be, ibe." (146.) 
5S, 5S, 458, 
eaT Gann e. Ta 


28. To illustrate by the same example the theory of gradually varying elements, let us 
establish the following definitions, for the present disturbed motion, 
k= Tl, ko=nNa— Wt, K,=N3— Wgt— wai 
A\y= 0}, Ag = We; A; = wz + gt, 
and let us call these six quantities «4, K2, K3, Ay, Ag, Ag the varying elements of that motion, by 
analogy to the six constant quantities e4, e3, 3, P1» Pe, P3, Which may, for the undisturbed 
motion, be represented in a similar way, namely, by (127.) and (128.), 
C= Tt, C:=Ng—Dat, €3=13—Wgt— dg, | 
P=, P2= Dp, Pz =0;+gt. 
We shall then have rigorously, for the six disturbed variables 7, , n2, N3, W1, D2, Wz, EXpressions 
of the same forms as in the integrals (127.) and (128.) of undisturbed motion, but with variable 
instead of constant elements, namely, the following: 
M=ki +A, Na=ka+ Aat, Ng =k +Azt— R 


(147.) 


(148.) 


(149.) 


and the rigorous determination of the six varying elements x, Ka, K3, A,, Àz, Àg, a8 functions of 
the time and of their own initial values e}, ¢,, €3, P1» Pe, Ps, depends on the integration of the 
6 following equations, in ordinary differentials of the first order, of the forms (105.): 


dk, oH, , 
d Whew» t (kı +A, 2), 
dk, SH, , 
dix, Bis.. ; 
a ere J 
and 
dà _ oH, x 
ie i u? (ky +A,t), 
a, ag ae 
ae Se = ae (Ka+ Aat), (151.) 
E oe: 


— — Á- a El 2 —— 2 
dt Skg á (kg +Agt zgi } 


H, being here the expression 


Hy =" ey AP (a HASH (cat Ast g, (152.) 


www.rcin.org.pl 


28, 29] II. GENERAL METHOD IN DYNAMICS 193 


which is obtained from (125.) by substituting for the disturbed coordinates 7, , na, na their values 
(149.), as functions of the varying elements and of the time. It is not difficult to integrate 
rigorously this system of equations (150.) and (151.); and we shall soon have occasion to state 
their complete and accurate integrals: but we shall continue for a while to treat these rigorous 
integrals as unknown, that we may take this opportunity to exemplify our general method of 
indefinite approximation, for all such dynamical questions, founded on the properties of the 
functions of elements C and E. Of these two functions either may be employed, and we shall 
use here the function C. 


29. This function, by (109.) and (152.), may rigorously be expressed as follows: 


pe t 
Cm | ote- ed+ age ddi 
ve (153.) 
+E [Ost ige- gjat: 


and has therefore the following for a first approximate value, obtained by treating the elements 
Ky, K2, Kg, Ay, Ag, Ag as constant and equal to their initial values e}, €z, €3, P15 Po» Ps; 


vat ara (ei + e3) + veji {p (pi + p3) + p53} 
a" "oh (154.) 
= 8 v“9pst+ 40 vg ‘ 
In like manner we have, as first approximations, of the kind expressed by the general formula 
(Z}.), the following results deduced from the equations (151.), 
A, =P1— p? (et+ 32,0), 
Az = P2— p? (Cat + $P2t), (155.) 
As=P—v* (ezt + 4p gt? — gt), 
and therefore, as approximations of the same kind, 


> 
e= —$p,t—- ae 


e= -ppt - anpe, (156.) 


Às— p 
e= -Pst + igt- 


Substituting these values for the initial constants e}, €, € in the approximate value (154.) for 
the function of elements C, we obtain the following approximate expression C’, for that func- 
tion, of the form supposed by our theory: 


_— _ 4 (A22)? + A- P2)? As- Ps)? 
Pe -g OA Seal opel | 
n {A1 — P1) Pa + (Àz — P2) Po + (Às — Ps) (Pa — 398) (157.) 
t? 2 2 2 2m2 t 2 & 272 
+ga H (Pit Pa) + D3} — ag” IPs + gorg. 


HMPII 25 


www.rcin.org.pl 


194 Ill. GENERAL METHOD IN DYNAMICS [29 


The rigorous function C must satisfy, in the present question, by the principles of the eighteenth 
number, the partial differential equation, 


8C_ pw? (/8C C 2) y2 /8C Eye 
T 5 (at +A t) + + (5, +t) }+ a p ist) (158.) 
and if it be put under the form (UŁ.), 
C = Ci + C, ’ 
C, being a first approximation, supposed to vanish with the time, then the correction C, must 
satisfy rigorously the condition 


_ a 8, , pw? (8C, 80, 50; o ji | ae 
C= {- ete Ott +t) + “(act tat) + =a +A,t—4gt 
1? òC, dC, 2 [8C2 
-afe (axe) +" (Sat) + (xe) J 
In passing to a second approximation we may neglect the second definite integral, and may 
calculate the first by the help of the approximate equations (155.); which give, in this manner, 


(159.) 


t 
0,= — | {O1—pr)*+ Qo pe)? + a-pa} dt 
p? fe 
+ p f {Ay (Ay — P1) + Àa (Àz — Pa) } dt 
y2 ft 
F z I. (As — 3gt) (Ag — pg) ° dt 
t 
nie {(Ay — Py)? + (Ag — P2)? + (Ag — P3)”} (160.) 
{8 > 
m 24 {u*p, (Ay — P1) + Be (Ag — P2) + v*Dz (Ag —Pps)} 
tt a DE Se Pale 
a Th g (\3— Ps) + zgo (# Pit Pps + vps 
ies Ay 4 + a 472 
240” IPs T 945" 9° 
We might improve this second approximation in like manner, by calculating a new definite 


integral C,, with the help of the following more approximate forms for the relations between the 
varying elements A,, Àz, A, and the initial constants, deduced by our general method: 


8C, C3 AP MoPi(14 we an -Pa (14 eee 


Sp, õp = pt 6" 24} 2 TI N A 


e = L001 80s UND E f h) 
a fop O prt Ss a) 12° 60)’ 


161. 
NE E E fa adm sk 
oe he Be AA OY BT ERY ID GO 
gè, TE Atty 
+F(14 tao) 
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in which we can only depend on the terms as far as the second order, but which acquire a correct- 
ness of the fourth order when cleared of the small divisors, and give then* 
Ay = Py — prt (e1 + 3P1t) + pË (e1 + 2744), 
Ag = Po— pt (e2 +$ Pat) + {eB (e2 + EPat), (162.) 
Àa = Pa — vt (e3 + $ Pat — ggt?) + ErP (e3 + Est — siogt?): 
But a little attention to the nature of this process shows that all the successive corrections to 
which it conducts can be only rational and integer and homogeneous functions, of the second 
dimension, of the quantities A,, Àa, As, P1; Pe, P3, 9, and that they may all be put under the 
following form, which is therefore the form of their sum, or of the whole sought function C; 
C= pa, (Ay— Py)? +5, 71 Ar— Pi) + WC, Pj 
+ pa, (Ag — Po)? +b, Do (Ag — Pa) + pC, D3 (163.) 
+ va, (As — ps)" + b, Ps (As — Ps) g ve, p3 ; 
+fr9 (Às — Pa) + hy gps + v’i,g?, 
the coefficients a, , &,, &c. being functions of the small quantities u, v, and also of the time, of 
which it remains to discover the forms. Denoting therefore their differentials, taken with 
respect to the time, as follows, 
da,=a,dt, da,=a,dt, &e., (164.) 
and substituting the expression (162.) in the rigorous partial differential equation (158.), we 
are conducted to the six following equations in ordinary differentials of the first order: 
2a, = (2a, + v*t)?; b, = (2a, + vt) (b, +t); = 3 (b, +t)’; | 
h= (2a, + vt) (f,— 3); h,= (b, +t) (= 31"); w=4(f,- $°)’; 
along with the 6 following conditions, to determine the 6 arbitrary constants introduced by 
integration, t 


(165.) 


ung stag AS aE E A A A 
te ioe me. | BO 
In this manner we find, without difficulty, observing that Ay» bi} c, may be formed from 
a,,b,,¢, by changing v to p, 
a,= —4v*t— $v cotan rt, a, = — pt — dp cotan pt, 


jis a (166.) 


l 1 t 
b,= —t+- „tan ni b = =f. — tan. 
2 p p 
1 pa | ut 
C=- atatan, aie = PT ba he 
i“ g (167.) 
f,= 28 — 4+ cotan vt, 
t2 t vt 
a n h 
7 42 
a oS tabi cotan vt. 


” Qy* 6y* 2,3 
* [If we substitute on the right-hand side of (150.) and (151.) ez, €z, €3, Pi» Pa» Pa for Ky, Kg, Kg, Az, Ag, Àg aS & 
first approximation, integration gives 
Ky =e +3 t (a +$pt), =p pt (e, +p), ete. 
The insertion of these values in the same equations followed by an integration as before gives (162.).] 
t [These coefficients come from (157.).] 


25-2 
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The form of the function C is therefore entirely known, and we have for this function of 
elements the following rigorous expression, 


Qai (Ay — 1)? + (Ag — P2)? _(As— Ps)? 
2u tan pt 2v tan vt 


a ; {A1 — Pi)? + (Ag— P2)? + (Ag — Pg)?} 
—t{p, (Ay — P1) + Po (Ag— P2) + P3 (Às — P3)} 


1 1 
ste Pata! (Ay — P1) + Po (Ay — P2)} tan Ká {as (Às — Pa) tan” (168.) 


' 
-z (i+ pi +p + (pi+p tan K + pitan” 


2 1 
+ i- m a “ cotan a) g (As — P3) + (5-; — -tan r) gps 


2 
+ (ss ee 5 cotan nt) g’, 


which may be variously transformed, and gives by our general method the following systems 
of rigorous integrals of the differential equations of varying elements, (150.), (151.): 


MOCT TDA TA _ Pry, 


KEE RET T rS 
i0. Py Pa at 

ey apy halt ps s tan z’ (169.) 
8C ng Àz — P3 fsi vt g t 1 

“3 Sp, vsinvt v is dre vt =) i 


50 1 1 
and So —(A,— pı) (14+ cotan pt) +p, (-1+/tan), 


1 
b0 1 1 pt 
== = —t+-— tan = 70. 
Ke Sa — (Ag —p,) (t+ © cotan ut) + pa ( Sa ant), (170.) 
80 1 1 vt se aa 
te e — (As — ps) (1+ <cotan wt) +p, ( -+5 tan 3) +9(5—at seotanst); 
that is, A, =p, cos pt —e, p Sin pt, l 
Ag = Pa COS pt — ep u SÌN pt, " (171.) 
Àz = Pz cos vt — ez vsin vt +g t- sinn), 
1 
and Ky = & (cos ut + utsin ut) +p, (7 sin pt — tos ut), 


j biz 
Ko = €p (COS ut + pet sin ut) + Po G sin ut — t cos mi) ; (172.) 


; Ei 2 
Ks = €z (cos vt + vt sin vt) + P3 (Fein vt — t cos v)-o (SS Z sinnt 5). 
y v 
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Accordingly, these rigorous expressions for the 6 varying elements, in the present dynamical 
question, agree with the results obtained by the ordinary methods of integration from the 6 
ordinary differential equations (150.) and (151.) and with those obtained by elimination from 
the equations (113.), (114.), (147.). 


Remarks on the foregoing Example. 
30. The example which has occupied us in the last six numbers is not altogether ideal, but 
is realised to some extent by the motion of a projectile in a void. For if we consider the earth as 


a sphere, of radius R, and suppose the accelerating force of gravity to vary inversely as the 
square of the distance r from its centre, and to be =g at the surface, this force will be repre- 


2 
sented generally by oo and to adapt the differential equations (78.) to the motion of a pro- 


jectile in a void, it will be sufficient to make 
v=o (2-5). (173.) 


If we place the origin of rectangular coordinates at the earth’s surface, and suppose the 
semiaxis of +z to be directed vertically upwards, we shall have 


r=V(R+2)?+2%+y?, (174.) 
and 
Hine 2 gee. Oar +) 
U + oR? (175.) 


neglecting only those very small terms which have the square of the earth’s radius for a divisor: 
neglecting therefore such terms, the force-function U in this question is of that form (110.) on 
which all the reasonings of the example have been founded; the small constants p, v being the 


real and imaginary quantities J Z, | =e, respectively. We may therefore apply the results 


of the recent numbers to the motions of projectiles in a void, by substituting these values for 
the constants, and altering, where necessary, trigonometrical to exponential functions. But 
besides the theoretical facility and the little practical importance of researches respecting such 
projectiles, the results would only be accurate as far as the first negative power (inclusive) of 
the earth’s radius, because the expression (110.) for the force-function U is only accurate so far; 
and therefore the rigorous and approximate investigations of the six preceding numbers, 
founded on that expression, are offered only as mathematical illustrations of a general method, 
extending to all problems of dynamics, at least to all those to which the law of living forces 
applies. 


Attracting Systems resumed: Differential Equations of internal or Relative 
Motion; Integration by the Principal Function. 


31. Returning now from this digression on the motion of a single point, to the more im- 
portant study of an attracting or repelling system, let us resume the differential equations (A.), 


which may be thus summed up: 
. dtòH =X (dy ôw — dw ðn); (A?.) 
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and in order to separate the absolute motion of the whole system in space from the motions of 
its points among themselves, let us choose the following marks of position: 


3 _2.mg _E.my Z .mz 


maita t said = 176. 
“ym? in? oe ( ) 


and 

Éi= ti- En, Ni=Yi— Yn: ki=zi— 2n; (177.) 
that is, the 3 rectangular coordinates of the centre of gravity of the system, referred to an 
origin fixed in space, and the 3n — 3 rectangular coordinates of the n — 1 masses My , Mg, ... Mp4; 
referred to the nth mass m,, as an internal and moveable origin, but to axes parallel to the 
former.* We then find, as in the former Essay, t 


T=} (0/2 + y;2-+2,2) Em i 
1 ; 78. 
+42, .m (E4240) — sa {(B, me) (B, m E mh] OT? 
the sign of summation £, referring to the first n — 1 masses only; and therefore, 
1 (/87\? (&7\?  (8T\? 1(/s7\? (sT\?_ (8T \? 
1 ~ 335m (e) +r) +z) nle) + (Gn) Gr) ) a 
sah cs ST ‘(2 ikea i (179.) 
tom {| F) Sy) 80’ . 
If then we put for abridgement, 
SLOP a Ss, ame’ 
+ eee =m ’ 
pig dy ov Bcd! 
Y, mõ =m ? (180.) 
pei: gpl. .mý' 
‘mse’ xm ’ 
we shall have the expression 
H=} (2, +y, +2?) Lm + 32, .m (x+y)? +z?) 
1 BA 
+H {Œ, -ma,)t+ (2, my’)? + (E, .ma’)9}—U, ay 
of which the variation is to be compared with the following form of (A?.), 
dtòH = (dx, dx) — dx; dx, + dy, dy, — dy} dy, +dz, dz, — dz} 8z,) a (C2) 
+2, .m (dE dx, — dx, SÉ + dy dy; — dy, n+ dl dz, — dz, 8%), : 


in order to form, by our general process, 6n differential equations of motion of the first order, 
between the 6n quantities z,,y,,2,,%,,Y,,2,,€ N, É, £, Ys, z, and the time ¢. In thus taking 
the variation of H, we are to remember that the force-function U depends only on the 3n — 3 
internal coordinates é, ņ, ¢, being of the form 


U = my (mifi + Maf: gt. t MnsSn—1) | 


; (D2.) 
+My Mofi a + MiM first. +My a Miya teens 


* (m=z, N etc.] 


t [P. 127, equation (69.).] 
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in which f; is a function of the distance of m; from m, , and f; ; is a function of the distance of m; 
from m,, such that their derived functions or first differential coefficients, taken with respect 
to the distances, express the laws of mutual repulsion, being negative in the case of attraction; 
and then we obtain, as we desired, two separate groups of equations, for the motion of the 
whole system of points in space, and for the motions of those points among themselves; namely, 
first, the group 

dz, =x, dt, d&s 0,| 


dy,=y,dt, dy)= (181.) 
dz,=z'dt, dz'=0, 
and secondly the group 
d=(z,+ >E ma, dt dp Bt a 
+ My, ig + , m 8& 
E E iB ; ,_ 180 (182.) 
dn=(yi+ E, .my,) at, dy, = 75 dt 
id i <3 war 
dč=(z; +E, .mz) dt d= g 


The six differential equations of the first order, (181.), between x, ,¥,,2,, x), Y, Z, and t, con- 
tain the law of rectilinear and uniform motion of the centre of gravity of the system ; and the 
6n — 6 equations of the same order, (182.), between the 6n — 6 variables £é, n, ¢, x) , y’ ,z’ and the 
time, are forms for the differential equations of internal or relative motion. We might eliminate 
the 3n — 3 auxiliary variables x’ . y’ , z, between these last equations, and so obtain the following 
other group of 3n — 3 equations of the second order, involving only the relative coordinates and 
the time, 


EM +4 00 
oman BE tm,” BE” 
wai, yU 
"aie tale: h (183.) 
("= 180 1 ôU. 
môt m, ' C’ 


but it is better for many purposes to retain them under the forms (182.), omitting, however, for 
simplicity, the lower accents of the auxiliary variables x’ , y’ , z; , because it is easy to prove that 
these auxiliary variables (180.) are the components of centrobaric velocity,* and because, in 
investigating the properties of internal or relative motion, we are at liberty to suppose that the 
centre of gravity of the system is fixed in space, at the origin of x,y,z. We may also, for simplicity, 
omit the lower accent of £,, understanding that the summations are to fall only on the first 
n— 1 masses, and denoting for greater distinctness the nth mass by a separate symbol M; and 
then we may comprise the differential equations of relative motion in the following simplified 


formula, ` 
dtòH =£ .m (dé òx’ Poly tia — dy’ 3n + dldz' — dz’ df), (E?.) 


in which H=}2z. m(z'24y2+z + yi: ma’)? + (X. my’)? + (X .mz’)*} — U. (F2.) 


* [By centrobaric velocity is meant the velocity relative to the mean centre: 
x, — x), = (1—2) + (a), —2/,) =f; Eme /Em=2;,.] 
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And the integrals of these equations of relative motion are contained (by our general 
method) in the formula 
8S =X. m (x'dE —a'da+ y’dn —b'5B +2'86 —c'dy), (G?.) 
in which g, B, y, a’, b’, c’ denote the initial values of £, n, f, x’, y’, 2’, and S is the principal function 
of relative motion of the system; that is, the former function S, simplified by the omission of the 
part which vanishes when the centre of gravity is fixed, and which gives in general the laws of 
motion of that centre, or the integrals of the equations (181.). 


Second Example: Case of a Ternary or Multiple System with one Predominant 
Mass; Equations of the undisturbed motions of the other masses about this, in 
their. several Binary Systems; Differentials of all their Elements, expressed by 
the coefficients of one Disturbing Function. 

32. Let us now suppose that the n — 1 masses m are smallin comparison with the nth mass M; 
and let us separate the expression (F?.) for H into the two following parts, 
H,== Dag) (a’2 +y'2+2z'?)— MX. mf, 

Mı Ma 
M 
Mi My 
M 
of which the latter is small in comparison with the former, and may be neglected in a first 
approximation. Suppressing it accordingly, we are conducted to the following 6n — 6 differential 
equations of the 1st order, belonging to a simpler motion, which may be called the undisturbed: 

ESAE dz’ 13H, _ yY. 


A, = (zita + YiY2 + %422— Mf, 2) +... | (H2.) 


a 


(Eir +Y Yn + ar — fin) +--+, 


di mse \ 'M d m Sé SE’ 
dy: 188, (,,m).,, dy _ 18H, N, : 
Bam by | +7)" dt m By S as 


dé 18H, m) , . de 18H, òf 
Eom (i+g) ii! Ba E Ol ae 
These equations arrange themselves in n—1 groups, corresponding to the »—1 binary 
systems (m, M); and it is easy to integrate the equations of each group separately. We may 
suppose, then, these integrals found, under the forms, 
K= x (t, $. N, C; x’, y’,2’), v= x(t, g, > 33 dahil | 
A=x (i, €,9, 6,2’, 9,2), T= E, 2,0, E (K?.) 
pos x” (t, g, UE Giz’, y’, z’), w= x® (t, ©; N, C, 2’, y’, z), 
the six quantities x, À, p, v, T, w being constant for the undisturbed motion of any one binary 
system; and therefore the six functions x, y®, x®, +, x, +, or K, À, u, v, T, w, being such as 
to satisfy identically the following equation, 
òk 8«6H, ôrõH, dkdH, 8x dH, 85H, Sx 5H 
Onm 51+ BE Ba’ Ba’ BE * Bn By’ By By 8 Be Be Be? 
with five other equations analogous, for the five other elements A, u, v, T, w, in any one binary 
system (m, M). 
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33. Returning now to the original multiple system, we may retain as definitions the equa- 
tions (K?.), but then we can no longer consider the elements x;, Àj, His Vis Ti, w; of the binary 
system (m;, M) as constant, because this system is now disturbed by the other masses my; how- 
ever, the 6n — 6 equations of disturbed relative motion, when put under the forms 

dé 8H, 5H, dz’ 6H, òH, 


a ae gt “eet sane. aaa 
mula 6H, 5H, dy’ dH, 5H, 2 
a e k alee a ae T Che) 
df 5H, 5H, dz’ òH, 5H, 
ap n “dlaehitel” ania ° Salat de 
and combined with the identical equations of the kind (L?.), give the following simple expression 
for the differential of the element «x, in its disturbed and variable state, 
os _o« dH, dk dH, d«dH, d« dH, kH, dx dH, N? 
=E Bc’ Sm BE tSn By’ — By’ By S Be Be! BL’ Sil 
together with some expressions for the differentials of the other elements. And if we 
express é, 1, %, x’, y’, z’, and therefore H, itself, as depending on the time and on these varying 
elements, we may transform the 6n — 6 differential equations of the 1st order, (M?.), between 
é, n, ¢,2’, y’, 2’, t, into the same number of equations of the same order between the varying 
elements and the time; which will be of the forms 


mo ={x, Nay +k, Hs, tte ak iA 
maA A EHA EHA PSR, to S 


mH ty, Wh iu, N Hiha N H ya, wi Xe 
dv òH 5H òH 8H òH ; 
m go Kh +h, Sith, a aa Pa teog 


dr ôH òH ôH 6H SH 
"k Sde; Kho +f, Nyt i ra -a + fr, wh, 


ôH ôH ôH 
mË fa, «} oe ? + {w AY y Heng tes tee 


if we put, for abridgement, 


ts B= 32 507 Ba’ 8E * By By! By/ By * BL Be’ Be’ 3L’ ra 
and form the other symbols {x, u}, {A,«}, &c., from this, by interchanging the letters. It is 
evident that these symbols have the properties, 

{A, K}= — {k, A}, {k, K}=0; (184.) 
and it results from the principles of the 15th number, that these combinations {«,A}, &c., when 
expressed as functions of the elements, do not contain the time explicitly. There are in general, 
by (184.), only 15 such distinct combinations for each of the n —1 binary systems; but there 
would thus be, in all, 15n—15, if they admitted of no further reductions: however, it results 


HMPII 26 
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from the principles of the 16th number, that 12n — 12 of these combinations may be made to 
vanish by a suitable choice of the elements. The following is another way of effecting as great a 
simplification, at least for that extensive class of cases in which the undisturbed distance 
between the two points of each binary system (m, M) admits of a minimum value. 


Simplification of the Differential Expressions by a suitable choice of the Elements. 
34. When the undisturbed distance r of m from M admits of such a minimum gq, corre- 
sponding to a time 7, and satisfying at that time the conditions 
=O, r”>0, (185.) 
then the integrals of the group (I?.), or the known rules of the undisturbed motion of m about 
M, may be presented in the following manner:* 


r=y{(éy' — a’)? + (ne! — Ly’)? + (tx' — éz)’; 


A=K— by’ +72’; 
M+m 19 12 12 ý 
u= zy (+y? +2) — Mf (r); 
ia en ne’ — Cy’ 
v=tan El — Cat”? | 
r M dr 
‘ENS 2 
a eat ae } (Q?.) 
E 2 +2Mf e) (1+) 5} 
ë 5 M) r? 
r M+m dr x 
—_— .—= .=.-dr 
kýr! M V dr? r 


w =v + sin t M u; 
V Bepa {+ 2Mf (r)— (1 +a st 

the minimum distance q being a function of the two elements x, u, which must satisfy the 
conditions way th m\ k? 
2u-+ 2Mf(q)— (1+ 37) a= Mf’ (a)+ (1+ 57) >; (186.) 

and sin-!s, tan-!t, being used (according to Sir John Herschel’s notation)t to express, not the 
cosecant and cotangent, but the inverse functions corresponding to sine and cosine, or the arcs 
which are more commonly called arc (sin=s), arc (tan =t). It must also be observed that the 


dr 
factor ——, 
ctor = 


is designed to be taken as equal to positive or negative unity, according as dr is positive or 
negative; that is, according as r is increasing or diminishing, so as to make the element under 
each integral sign constantly positive. In general, it appears to be a useful rule, though not 


always followed by analysts, to employ the real radical symbol V R only for positive quantities, 
r 


Vr 
* [See Appendix, Note 4, p. 623, and also footnote on p. 206.] 
+ [For the history of the notation for the inverse trigonometrical functions see Cajori, History of Mathematical 
Notations, Vol. 1, pp. 175-178.] 


which we have introduced under the signs of integration, is not superfluous, but 


unless the negative sign be expressly prefixed; and then will denote positive or negative 
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unity, according as r is positive or negative. The arc given by its sine, in the expression of the 
element w, is supposed to be so chosen as to increase continually with the time. 

35. After these remarks on the notation, let us apply the formula (P?.) to calculate the 
values of the 15 combinations such as {x, À}, of the 6 constants or elements (Q?.). 

Since 


raV(e+nt+e), (187.) 
it is easy to perceive that the six combinations of the 4 first elements are as follows: 
{x,A}=0, {x,p}=0, {x,vy}=0, fA,pw}=0, fA,y}=1, {y,v}=0. (188.) 


To form the 4 combinations of these 4 first elements with 7, we may observe, that this 5th 
element 7, as expressed in (Q?.), involves explicitly (besides the time) the distance r, and the 
two elements «, p; but the combinations already determined show that these two elements may 
be treated as constant in forming the four combinations now sought; we need only attend, 
therefore, to the variation of r, and if we interpret by the rule (P*.) the symbols {«x, r}, {A,r}, 
{m.r}, {v,r}, and attend to the equations (I?.), we see that 


{x,r}=0, {A,r}=0, rE NR a {v,r}=0, (189.) 


A being the total differential coefficient of r in the undisturbed motion, as determined by the 
equations (I?.); and, therefore, that 

{x,7}=0, fA,7}=0, {v,7}=0, (190.) 
and 


{u,7}= he oe ™ i (191.) 


observing that in differentiating the expressions of the elements (Q?.), we may treat those 
elements as constant, if we change the differentials of £, 7, £, x’, y’, z’ to their undisturbed values. 
It remains to calculate the 5 combinations of these 5 elements with the last element w; which is 
given by (Q?.) as a function of the distance r, the coordinate ¢, and the 4 elements x, A, p, v; so 
that we may employ this formula, 


dw Sw 8 dw dw dw 
{e, w= {e, h+ {e, +z ie k} + SÀ {e, À} + Sp fe, H+ {e, v}, (192.) 


in which, if e be any of the first five elements, or the distance r, 


ô ô 
{e, r}= 1g tage Hos, S) {e, {j= -— {e,x}=0, (193.) 


Sw (ðk)! Sw dew ðw .. 
se-la) gaa aijo 


the formula (192.) may therefore be thus written: 


and 


(195.) 


,(. Se de Se 
rp (Sse +59 Sse) sel faoa 
BAY rae ie] (se) 
He} 3 GM +5" leh 


26-2 


www.rcin.org.pl 


204 HMHI. GENERAL METHOD IN DYNAMICS [35 


We easily find, by this formula, that 
dr òw. 


{x,w}=—-1; fA,a}=0; {u,w}=0; f{r,w}= dt Sp’ (196.) 
and 
dv ðw dw 
Yo} = — 55 er — gy 0 (197.) 


The formula (195.) extends to the combination {7,w} also; but in calculating this last 
combination we are to remember that 7 is given by (Q?.) as a function of x, p, r, such that, 


br dt 


and thus we see, with the help of the combinations (196.) already determined, that 
òr w 98 [* 68 ad 
ag EEE MEY: Ordr; | 2dr, — 


if we represent for abridgement by ©, and Q, the coefficients of dr under the integral signs in 
(Q?.), namely, 


M = m K?\-Ż 
rer i shag (200.) 
M+m ar F adh K 
These coefficients are evidently connected bat the relation 
| 80,  5Q, 
Set 37 (201.) 
which gives 
8 a 6? 
se|, tr, |, nro, Hien 


r, being any quantity which does not vary with the elements «x and u; we might therefore at 
once conclude by (199.) that the combination {7, w} vanishes, if a difficulty were not occasioned 
by the necessity of varying the lower limit q, which depends on those two elements, and by the 
circumstance that at this lower limit the coefficients ©,, Q, become infinite. However, the 
relation (202.) shows that we may express this combination {r, w} as follows: 


5 fr 5 f7 
{r,0}=2. j @,dr+ 5 Í “yar, (203.) 


r, being an auxiliary and arbitrary quantity, which cannot really affect the result, but may be 
made to facilitate the calculation; or in other words, we may assign to the distance r any 
arbitrary value, not varying for infinitesimal variations of x, u, which may assist in calculating 
the value of the expression (199.). We may therefore suppose that the increase of distance r — q 
is small, and corresponds to a small positive interval of time t—7, during which the distance r 
and its differential coefficient 7’ are constantly increasing; and then after the first moment 7, 
the quantity 


l 
0,=7 hi (204.) 
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will be constantly finite, positive, and decreasing, during the same interval, so that its integral 
must be greater than if it had constantly its final value; that is, 


t-r= |" 0,ar> (r—q),. (205.) 
a 


Hence, although ©, tends to infinity, yet (r —q)©, tends to zero, when by diminishing the 
interval we make r tend to q; and therefore the following difference 


foa AAN @,dr -1T (5-5) 9dr (206.) 


will also tend to 0, and so will also its partial differential coefficient of the first order, taken with 
respect to u.* We find therefore the following formula for {r, w}, (remembering that this com- 
bination has been shown to be independent of r,) 


froa}= A {> -Soari TTR = "|. earls (207.) 


the sign.. A _ implying that the limit is to be taken to which the expression tends when r tends 


to q. Ea as in (199.), the integral |” @,dr may be considered as a known 


function ofr, q, x, p, or simply ofr, q, x, if p be eliminated by the first condition (186.); and since 
it vanishes independently of x when r =q, it may be thus denoted: 


['ear— (r, q, K) ime $ (q, q, k), (208.) 


the form of the function ¢ depending on the law of attraction or repulsion. This integral there- 
fore, when considered as depending on x and p, by depending on « and q, need not be varied with 
respect to x, in calculating {7, w} by (207.), because its partial differential coefficient l f @, dr) ; 
q 
obtained by treating q as constant, vanishes at the limit r =q; nor need it be varied with respect 
to q, because, by (186.), 
ôg _Mim«sq_,. 
a ago: (209.) 


it may therefore be treated as constant, and we find at last 
{r, w} =0, (210.) 


the two terms (199.) or (203.) both ere to infinity when r tends to q, but always destroying 
each other. 

* | Fir is of the form —¥ where ys remains finite, . -(4- = 5) is of the form . 
the interval of integration. It follows that a (4-4) 0 „dr > 0 as r>a | 


J where fis finite in 
/r—q 


www.rcin.org.pl 


206 II. GENERAL METHOD IN DYNAMICS [36 


36. Collecting now our results, and presenting for greater clearness each combination 
under the two forms in which it occurs when the order of the elements is changed, we have, for 
each binary system, the following thirty expressions:* 

{x, }=0, {x, E {x, v}=0, {x,7}=0, {x,w}=—-1, 

fA, K}=0, fA, p}=0, {A, vy}=1, {A, 7}=0, fA, w}=0, 

{u, }=0, {w,A}=0,  {u, =0,  {u,7}=1, {u, w}=0, (R2.) 
{v, c}=0, {v, AJ=—1, {v, p}=0, {v, T}=0, {v, w}=0, 

{r, e}=0, {7, A}=0, {r, p}=—-1, {7, y}=0, {7, w}=0, 

{w,k}=1, {w,A}=0, {w, p} = 0, {w,v}=0, {w, T}=0; 

so that the three combinations 
{HT}, {a}, {A,r} 


are each equal to positive unity; the three inverse combinations 
TH} {Kw}, {vA} 
are each equal to negative unity; and all the others vanish. The six differential equations of the 
first order, for the 6 varying elements of any one binary system (m, M), are therefore, by (O2.), 
m- òH, dr ôH, 
git? ae ar E 
mee = dH, dx òH 
ar te? me io i adi 
mI 3H 4,47 3H, 
di Op. dt dA’ 
and, if we still omit the variation of t, they may all be summed up in this form for the variation 


of H,, 
8H,=X.m (p'òr— rT Òu +w Òr — k'8w + Av — v'à), (T?.) 
which single formula enables us to derive all the 6n — 6 differential equations of the first order, 
for all the varying elements of all the binary systems, from the variation or from the partial 
differential coefficients of a single quantity H,, expressed as a function of those elements. 

If we choose to introduce into the expression (T?.), for 5H,, the variation of the time t, we 
have only to change ôr to 5r—58t, because, by (Q?.), ôt enters a so accompanied; that is, t 
enters only under the form t—7;, in the expressions of é; ; 7;, bis Xi» Yi, 2; a8 functions o the time 
and of the elements; we have, therefore, 

ao y a Dm; (211.) 

* [If we put ¢=r cos l cos ¢, n=r cos l sin ¢, =r sin l, the principal function S satisfies 

2 2 2 
at a (æ) +(Fa) + Crona) J7 =o 


The Jacobi Complete Integral can be put in the form 


S=-ui+(-\) 4+] [a ay J ae [foros (145) SY art const 


where «x, A, p are arbitrary constants. We, then, easily obtain Ba, = —v; 08 —7r. Hence by Donkin’s 


Theorem (Routh’s Advanced Rigid Dynamics, Art. 496), (x, — w), (A, v), (u, 7) form a canonical set of constants and 
hence relations (R?.) hold.] 
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and since, by (H?.), (Q?.), 


H,=2.mp, (212.) 
we find finally, dH, 3H, 
netat e: 2 
/ ian i adk 


This remarkable form for the differential of H, , considered as a varying element, is general 
for all problems of dynamics. It may be deduced by the general method from the formulæ of 
the 13th and 14th numbers, which give* 


dH, dH. (8H, dx, _ dH, w) òH, 5 (SH Sken _ dH, een) 
dt BA 8 sa ôw òn ey Oa ôn 50 òw òn (213.) 
_ôH, 8x, , 5H, dK, SH Sken SH 


ie a c aia. at ee St” 


K1, Kg, ++» Kg, being any 6n elements of a system expressed as functions of the time and of the 
quantities 7, w; or more concisely by this special consideration, that H, + H, is constant in the 
disturbed motion, and that in taking the first total differential coefficient of H, with respect to 
the time, the elements may by (F'.) be treated as constant. It is also a remarkable corollary of 
the general principles just referred to, but one not difficult to verify, that the first partial 
differential coefficient a of any element «,, taken with respect to the time, may be expressed 
as a function of the elements alone, not involving the time explicitly. 


On the essential distinction between the Systems of Varying Elements considered in 
this Essay and those hitherto employed by mathematicians. 


37. When we shall have integrated the differential equations of varying elements (S*.), we 
can then calculate the varying relative coordinates é, 7, ¢, for any binary system (m, M), by the 
rules of undisturbed motion, as expressed by the equations (I*.), (Q?.), or by the following 
connected formule: ' 


=r (cos 0+ 2sin (0—»)siny), 

=r sin 0 LS TA ) cos ) 2 
g ies k v)COSV], (V .) 
t= V2 M sin (0—7): 


* [To deduce (213.) from §§ 13, 14. We have (since H, remains constant for undisturbed motion) 


dH, _» OH, Ong dri See 2G ae ta! ae) 8 oH, 
dt On, Ox; dt 0m, Ox; dt Ons Ox; OK 4 Ox; 

Sadap a (Oe d d et) GE, 

Ons x; tom, Ox; One Om; Om; Ont OK; 

-yH oH, {71 oe Ox; Ox; _ Ons OK; ad 4s oH, fone s OK; Ox; _ oD, OK; oh |. 
Ox; Ons OK; Ont am, ~ Ox; Om, Ont an, OK; One Om, OK; 0m; On 
Ons One _ Ons Ones _ OW, Ox; Om, Ox; 
But & Bx; ng O a bua, 0, Dki On =0, Dx; Fd, =5,,, where ô equals 1 or 0 according as s is or is 
not equal to te 3 dH, _ y OH; oH, care On; _ 0H, ett oy y Oe OK; J 
aed OK; On; OW; Om, Ons OK, ôt 
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in which the distance r is determined as a function of the time ¢ and of the elements 7, «, p, by 
the 5th equation (Q?.), and in which* 


0=w+ (W?.) 


2)? 
Jermo- ga 
q being still the minimum of r, when the orbit is treated as constant, and being still connected 
with the elements «x, p, by the first equation of condition (186.). In astronomical language, M 
is the sun, m a planet, £, ņ, { are the heliocentric rectangular coordinates, r is the radius vector, 
0 the longitude in the orbit, w the longitude of the perihelion, v of the node, @—w is the true 
anomaly, —v the argument of latitude, p the constant part of the half square of undisturbed 
heliocentric velocity, diminished in the ratio of the sun’s mass (M) to the sum (M +m) of masses 


of sun and planet, « is the double of the areal velocity diminished in the same ratio, a is the 


versed sine of the inclination of the orbit, q the perihelion distance, and 7 the time of perihelion 
passage. The law of attraction or repulsion is here left undetermined; for Newton’s law, p is 
the sun’s mass divided by the axis major of the orbit taken negatively, and « is the square root 
of the semiparameter, multiplied by the sun's mass, and divided by the square root of the sum 
of the masses of sun and planet. But the varying ellipse or other orbit, which the foregoing 
formule require, differs essentially (though little) from that hitherto employed by astronomers: 
because it gives correctly the heliocentric coordinates, but not the heliocentric components of 
velocity, without differentiating the elements in the calculation; and therefore does not touch 
but cuts, (though under a very small angle,) the actual heliocentric orbit, described under the 
influence of all the disturbing forces. 


38. For it results from the foregoing theory, that if we differentiate the expressions (V?.) 
for the heliocentric coordinates, without differentiating the elements, and then assign to those 
new varying elements their values as functions of the time, obtained from the equations (S?.), 
and deduce the centrobaric components of velocity by the formule (I?.), or by the following: 

Tg ,_. M7’ 7 ae 
, M =a y TET *=M+m’ (3144 
* [To deduce (V?.) from (Q?.) we have 
(nz — Ly’)? +(x — Ez) =x? —(«—A)? =(2xA—D?). 


Therefore ne’ — by’ =V2eA—M sin v, Ez’ — Lx’ =V 2eA—? cos v. 

Consequently V 2xA—® (£ sin v—n cos v) +(«—A) [=LE (nz —Ly’)=0, 

that is, g sin yn cos y= — zt. 

Also Etere. | 
Hence, Ẹ cos v +y sin v-e- of (e-a 

From (W?.) and (Q*.) we see that 7 =r sin (@—v), and therefore 


é sin y—7 cos y=(-142) rsin (@—v), 


€ cos v +7 sin v=r cos (0—v), 
from which (V*.) follows immediately.] 
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then these centrobaric components will be the same functions of the time and of the new varying 
elements which might be otherwise deduced by elimination from the integrals (Q?.), and will 
represent rigorously (by the extension given in the theory to those last-mentioned integrals) the 
components of velocity of the disturbed planet m, relatively to the centre of gravity of the whole 
solar system. We chose, as more suitable to the general course of our method, that these centro- 
baric components of velocity should be the auxiliary variables to be combined with the helio- 
centric coordinates, and to have their disturbed values rigorously expressed by the formule 
of undisturbed motion; but in making this choice it became necessary to modify these latter 
formule, and to determine a varying orbit essentially distinct in theory (though little differing 
in practice) from that conceived so beautifully by Lagrange.* The orbit which he imagined was 
more simply connected with the heliocentric motion of a single planet, since it gave, for such 
heliocentric motion, the velocity as well as the position; the orbit which we have chosen is 
perhaps more closely combined with the conception of a multiple system, moving about its 
common centre of gravity, and influenced in every part by the actions of all the rest. Whichever 
orbit shall be hereafter adopted by astronomers, they will remember that both are equally fit 
to represent the celestial appearances, if the numeric elements of either set be suitably deter- 
mined by observation, and the elements of the other set of orbits be deduced from these by 
calculation. Meantime mathematicians will judge, whether in sacrificing a part,of the simplicity 
of that geometrical conception on which the theories of Lagrange and Poisson are founded, a 
simplicity of another kind has not been introduced, which was wanting in those admirable 
theories; by our having succeeded in expressing rigorously the differentials of all our own new 
varying elements through the coefficients of a single function: whereas it has seemed necessary 
hitherto to employ one function for the Earth disturbed by Venus, and another function for 
Venus disturbed by the Earth. 


Integration of the Simplified Equations, which determine the new varying Elements. 


39. The simplified differential equations of varying elements, (S?.), are of the same form as 
the equations (A.), and may be integrated in a similar manner. If we put, for abridgement, 
(r, e= aee v t) Ha} dt, (X2.) 

and interpret similarly the symbols (u, w, A), &c., we can easily assign the variations of the 

following 8 combinations, (T, K, v), (H, w, A), (u, K, v), (T, w, A), (T, w, v), (m, £, À), (T, K, À), (u, w, v); 

namely, 

ò (T, K, v) =E .m (Thu — To Spry + KÕw — Ky dwg + v SÀ — vo dào) — Hy Òt, 

8 (u, w, A) =X. m (U9 O79 — LÒT + wg Sky — WSK + Ap dr — Adv) — H, dt, 

8 (u, K, v) =X. mM (wy S79 — Sr + Kdw — kowo +v ÖA — vodà) — H, dE, 

8 (T, w, A) =z . m (7 dpe — To Spy + wo Öko — WSK + Ag drvg — Adv) — H, dt, 


8 (T, w, v) =X. m (Thu — Ty Spry + Wy Sky — WSK + VbA— vo ðo) — H zòt, ity 
Ò (u, K, A) =X. mM (pg S79 — HÒT + KOwW — Ky dw + Ay dv — Adv) — H, dt, 
ô (T, K, A)=U.m (Tõu — Ty) Sprg + K Bw — Ky Sw + Ay Svq — Adv) — H, dE, 
8 (u, w, v) =È . M (wg S79 — OT + Wy Sky — WSK + VSA— vodà) — H zòt, 
| * [See Appendix, Note 7, p. 628.] 
HMPII 27 
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Ko; Ao» Hos Yos To» wo being the initial values of the varying elements «x, A, u, v, T, w. If, then, we 
consider, for example, the first of these 8 combinations (r, x, v), as a function of all the 3n — 3 
elements u;, w;,A;, and of their initial values jp ;, w9,;, Ao,;, involving also in general the time 
explicitly, we shall have the following forms for the 6n—6 rigorous integrals of the 6n—6 
equations (S?.): 


ô ò 
mirim OE) Miroi — 5 (yk, | 


5 8 

Mie Da (tHe); MiKo, i = pE PEA e. K; hae (Z?.) 
5 8 

hs TT (T, Kv); Mivo, i = ars TTH (T, K, v); 


and in like manner we can deduce forms for the same rigorous integrals, from any one of the 
eight combinations (Y*.). The determination of all the varying elements would therefore be 
fully accomplished, if we could find the complete expression for any one of these 8 combinations. 


40. A first approximate expression for any one of them can be found from the form under 
which we have supposed H, to be put, Paa as a function of the elements and of the time, 
which may be thus denoted: 


H,=H, (t, K1, Àr Ha: Vis Tis W15 +e. Kn—1> E h E Vn—1>Tn-1> wn); (A3.) 


by changing in this function the varying elements to their initial values, and employing the 
following approximate integrals of the equations (S?.), . 


rine pdt fide vat adh 8H, 
wpa += | Pat kiia | cana —— dt, 


0 m J o Spo 
1 f'dH, 1 [73H 7 
Om Pete en, E Oy ae ra (B .) 
j, 1 (ôH, 1 f'dH, | 
A rote be dt, yan — | Scat 
For if we denote, for example, the first of the 8 combinations (Y?.) by G, so that 
G=(r,x,v), (o) 
we shall have, as a first approximate value, 
t ôH, òH, òH, 
G,= {2 g Sn, E m Sr, 24% 5 ) H, dt; (D?.) 


and after thus expressing G} as a function of the time, and of the initial elements, we can 
eliminate the initial quantities of the forms 79, Ko, vo, and introduce in their stead the final 
quantities u, w, A, so as to obtain an expression for G4 of the kind supposed in (Z?.), namely, a 
function of the time t, the varying elements p, w, A, and their initial values pip , wg, Ag. An approxi- 
mate expression thus found may be corrected by a process of that kind, which has often been 
employed in this Essay for other similar purposes. For the function @, or the combination 
(r, x, v), must satisfy rigorously, by (Y?.), (A®.), the following partial differential equation: 


dG 1 òG 1864 184 1 8G ) 
Wn- 


0=— H t — —,; , ARA i a E 4S ea. PIE bi, y 
T ( Mı bw, . Pa Mı dA, m, Spy 1 Me Swa - ) 


òt 
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and each of the other analogous functions or combinations (Y?.) must satisfy an analogous 
equation: if then we change G to Gh + Gh, and neglect the squares and products of the coeffi- 
cients of the small correction G,, G} being a first approximation such as that already found, we 
are conducted, as a second approximation on principles already explained, to the following 
expression for this correction G: 
GD td 1 8G, 1 8G, 1 84, $ % 
G,= -fi ae +H, (. m, dw,’ A> His m, dA , m, õm’ Wi; .-) fat: (F .) 

which may be continually and indefinitely improved by a repetition of the same process of 
correction. We may therefore, thecretically, consider the problem as solved; but it must remain 
for future consideration, and perhaps for actual trial, to determine which of all these various 
processes of successive and indefinite approximation, deduced in the present Essay and in the 
former, as corollaries of one general Method, and as consequences of one central Idea, is best 
adapted for numeric application, and for the mathematical study of phenomena.* 


* [Two memoirs by M. Houél presented to the Faculté des Sciences de Paris (1855) are of interest. In the first 
he employs the Principal Function to prove various known theorems about planetary perturbations and in the 
second applies this method to detailed calculations of the perturbations of Jupiter. A complete bibliography of 
works relating to Hamilton’s dynamical methods prior to that of Houél was given by Cayley, Brit. Ass. Report 
(1857), p. 40.] 
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